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I  Introduction 


l.l  Objectives 

In  theory,  a  Space-Time  Adaptive  Processor  (STAR)  based  radar  system  can  be  used  to  simulta¬ 
neously  suppress  clutter  and  wideband  noise  jammers  (WNJ),  Simultaneous  suppression  typically  involves 
the  estimation  of  a  joint  space-time  covariance  matrix  which  can  be  quite  large  even  for  moderately  sized 
radar  systems  required  to  handle  a  number  of  WNJ.  Since  the  statistics  of  the  interference  sources  generally 
must  be  estimated,  this  approach  may  be  suboptimal  in  terms  of  both  performance  and  cost.  The  size  of  the 
secondary  data  set  required  to  adequately  estimate  the  covariance  matrix  may  be  unattainable,  particularly 
in  a  nonhomogenous  interference  environment.  Additionally,  since  WNJ  have  tem[X)raIIy  white  speclrums, 
Doppler  processing  of  WNJ  is  ineffective  and  causes  an  unnecessary  inflation  of  the  size  of  the  space-time 
covariance  matrix  to  be  estimated.  An  approach  to  the  problem  of  handling  combined  clutter  and  WNJ 
interference  is  to  panition  the  adaptive  processor  so  that  clutter  and  WNJ  can  be  handled  separately.  A 
promising  approach  already  addressed  in  the  literature,  which  provides  the  desired  partition,  makes  use  of 
the  concept  of  jammer  pre-suppression  followed  by  a  clutter  rejection  filter. 

This  report  addresses  the  topic  of  Spatial  Adaptive  Pre-Suppression  (SAPS)  of  WNJ  in  conjunction 
with  STAR  Specifically,  system  architectures  which  cascade  an  Adaptive  Spatial  Pre-Processor  (ASP)  for 
WNJ  suppression  with  STAR  are  investigated.  The  objective  of  this  study  is  to  evaluate  these  cascaded  sys¬ 
tems  as  alternatives  to  simultaneous  suppression,  for  the  case  when  the  available  secondary  data  set  is  lim¬ 
ited,  and  to  elucidate  some  of  the  key  issues  in  determining  the  most  suitable  approach. 

1.2  Overview 

This  report  extends  the  original  work  described  by  Klemm  in  (2J,  and  later  on  from  a  slightly  dif¬ 
ferent  perspective  by  Marshall  in  [3J,  pertaining  to  the  cascaded  (or  two-step  nulling)  pre-suppression  of 
WNJ’s  in  STAR  based  radar  systems.  Two  aspects  of  the  cascaded  ASP  /  STAR  approach  are  addressed:  ( I ) 
system  architecture  and  it’s  impact  on  detection  performance,  and  (2)  ASP  algorithm  development  and  the 
acquisition  of  clutter-free  jammer  secondary  data. 

The  underlying  theme  behind  the  SAPS  approach  (i.e.  reduction  in  the  size  of  the  statistical  cstima- 
tion  problem)  is  also  one  motivation  behind  the  partially  adaptive  techniques  which  have  been  under  inves- 
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tigaiion  since  the  I970’s  in  the  case  of  partially  adaptive  arrays  and  since  the  late  I980's  in  the  case  of 
partially  adaptive  STAP.  The  benefits  of  degrec-of- freedom  (DOF)  reduction  in  terms  of  detection  perfor¬ 
mance  and/or  system  affordability  can  be  substantial.  As  a  result,  the  cascade  architectures  proposed  herein 
are  a  combination  of  partially  adaptive  array  techniques  and  a  high  performance  partially  adaptive  STAP 
approach.  The  idea  of  separating  the  WNJ  and  clutter  suppression  problems  provides  a  new  perspective  in 
the  general  area  of  partially  adaptive  processing.  By  taking  advantage  of  the  inherent  differences  between 
the  WNJ  and  clutter  interference,  a  large  covariance  estimation  problem  is  reduced  to  two  much  smaller 
ones.  In  Section  2,  three  cascaded  architectures  and  a  baseline  ASP  algorithm  are  proposed  and  compared 
with  two  simultaneous  STAP  approaches.  Throughout  this  report,  computer  generated  analytical  results  are 
used  as  a  means  for  performance  evaluation  and  comparison.  Therefore,  part  of  Section  2  is  devoted  to  a 
description  of  the  computer  models  used  as  well  as  the  performance  metrics.  Section  2  results  are  generated 
assuming  the  availability  of  a  sufficient  number  of  jammer-only  secondary  data  samples  to  accurately  esti¬ 
mate  the  Jammer  covariance  matrix.  The  problem  of  obtaining  these  secondaries  is  the  subject  of  Section  3. 

A  key  issue  in  applying  the  SAPS  concept  is  the  acquisition  of  the  necessary  jammer-only  statistics 
for  adaptive  suppression,  free  from  strong  clutter  contamination.  Common  acquisition  methods  include  the 
use  of  clutter-free  range  gates  for  low  pulse-repetition-frequency  (PRF)  systems,  clutter-free  Doppler  bins 
for  high  PRF  systems,  or  receive-only  mode  between  two  coherent  processing  intervals  (CPI).  All  of  these 
techniques  require  jammer  data  to  be  collected  within  a  restricted  region  of  the  available  space-time 
domain,  and  may  not  always  be  able  to  generate  sufficient  jammer-only  data.  Moreover,  fast-changing  jam¬ 
ming  environments  and  large-scale  PRF  hopping  can  also  make  these  techniques  unsuitable.  In  Section  3, 
a  new  technique  is  presented  which  makes  use  of  a  carrier  frequency  sideband,  close  to  but  disjoint  from 
the  radar's  mainband,  to  estimate  the  jammer-only  covariance  matrix.  Such  an  idea  can  be  applied  to  a  sys¬ 
tem  with  any  PRF,  and  the  entire  or  any  appropriate  portion  of  the  Range  Processing  Interval  (RPI)  could 
be  used  to  collect  jammer  data.  It  should  be  noted  that  wideband  jammers  are  designed  to  sufficiently  cover 
the  radar’s  mainband,  making  a  sideband  containing  their  energy  always  available.  Since  the  ASP  weights 
for  jammer  suppression  are  derived  from  data  with  a  carrier  frequency  different  from  the  radar  signal,  an 
important  part  of  the  successful  implementation  of  the  sideband  approach  involves  the  design  of  sideband 
compensation  algorithms.  These  and  other  algorithms  are  investigated  in  Section  3. 


2  Cascaded  Space  and  Space-Time  Adaptive  Processing 

TTie  intent  of  the  Cascaded  Space  and  Space-Time  Adaptive  Processor  (CS_STAP)  is  to  separate 
the  WNJ  suppression  function  from  the  clutter  suppression  function,  thereby  replacing  a  large  space-time 
covariance  matrix  estimation  problem  with  the  estimation  of  two  smaller  covariance  matrices.  Some  of  the 
potential  benefits  of  this  approach,  relative  to  a  system  employing  STAP  to  simultaneously  reject  jammers 
and  clutter,  are  (1)  improved  detection  performance  in  an  environment  of  WNJ's  and  nonhomogenous  clut¬ 
ter.  (2)  reduced  susceptibility  to  false  detections  by  discretes  in  the  same  range  cell,  and  (3)  a  more  cost- 
efficient/practical  system  architecture. 

2.1  Fundamental  System  Models 

In  the  following,  a  system  model  is  established  which  is  commonly  used  for  airborne  surveillance 
pulsed-Doppler  radars.  For  a  well  organized  description  of  a  pulsed-Doppler  radar  system  model  see  for 
example  [91.  An  abbreviated  development  of  the  models  peninent  to  this  study  are  given  in  Section  2.1.1 
thru  Section  2. 1 .6.  The  data  models  of  interest  here  will  be  complex  sampled  data  vector  descriptions  of  the 

received  radar  signals  along  with  their  I  **  and  2"^  order  statistics.  This  will  be  sufficient  to  establish  the  key 
performance  metrics  (e.g.  probability  of  detection,  Pj,  at  a  specified  probability  of  false  alarm,  Pf)  for  the 
case  when  the  data  is  Gaussian  distributed. 

2.1.1  Radar  System 

The  radar  system  considered  here  employs  a  phased-anay  antenna  which  provides  Ng  column  sen¬ 
sor  outputs  whose  phase  centers  correspond  to  a  uniformly  spaced  line  array  along  the  y-axis.  The  coordi¬ 
nate  system  used  is  depicted  in  Figure  2-1 .  A  unit  vector,  r,  pointing  in  the  (4>,6)  direction  is  given  by 


r  =  cos4»cos0-x-t- sin4>cos0-y +  sin6-z  (2-1) 

where  x.y,  and  z  are  unit  vectors.  Vectors  corresponding  to  physical  quantities  (i.e.  direction,  velocity,  etc.) 
will  be  distinguished  from  data  vectors  by  using  italics.  The  array  is  situated  so  that  the  first  element  resides 
at  the  origin  of  the  x-y  plane  and  along  the  z-axis  at  the  platform  altitude.  The  first  element  serves  as  the 
phase  reference  for  the  entire  array.  The  x-y  plane  is  tangent  to  the  earth’s  surface  directly  below  the  array. 
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} 


The  n**'  sensor  position  is  described  by  it’s  position  vector.  relative  to  the  first  sensor.  For  the  case  con 
sidered  here  in  Figure  2-1  we  have 


=  Rj  •  d  y  .  (2-2) 

The  angle  directions  are  defined  as  follows:  (1 )  positive  azimuth,  <)),  is  a  counter-clockwise  rotation  from 
the  x-axis  to  the  y-axis,  (2)  positive  elevation,  0.  is  from  the  horizontal  xy-plane  down  towards  the  negative 
z-axis.  and  (3)  positive  platform  velocity  azimuth  misalignment  angle,  4>p.  is  a  clockwise  rotation  from  the 
y-axis  to  the  x-axis.  The  platform  velocity  azimuth  misalignment  angle,  4>p,  is  the  angle  from  the  y-axis  to 
the  projection  of  the  platform  velocity  vector  onto  the  xy-plane.  The  numerical  results  generated  in  this 
report  will  always  have  the  platform  velocity  vector  residing  in  the  xy-plane  as  indicated  in  Figure  2-1.  so 
that  only  (fip  is  necessary  to  define  the  velocity  direction.  Although  the  results  presented  are  for  horizontal 
platform  motion,  the  formulas  involving  the  platform  velocity,  Vp,  are  for  a  general  vector.  The  line  array  is 
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of  great  practical  interest  in  radar  engineering  and  will  be  sufficient  in  this  investigation  to  establish  the 
desired  results.  The  reader  interested  in  a  detailed  treatment  of  an  array  with  arbitrary  sensor  positions  is 
referred  to  (10). 

A  top-level  block  diagram  of  the  radar  system  considered  is  given  in  Figure  2-2.  For  convenience. 


Rgure  2-2:  Radar  Signal  Processor  Top-Level  Block  Diagram 


the  receiver  (i.e.  down-conversion  to  baseband,  waveform  matched  filter,  and  sampling  and  digitization  via 
A/D  conversion)  portion  of  the  radar  system  is  imbedded  in  the  Adaptive  Spatial  Processor  (ASP)  block. 
Different  spatial  processing  architectures  can  have  different  receiver  schemes  in  order  to  minimize  system 
cost  and  control  channel  errors.  Therefore,  the  receiver  structure  for  each  will  be  described  separately  as 
needed.  Narrowband  signals  are  considered  in  this  paper,  unless  stated  otherwise,  and  therefore  the  analog 
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signals  prior  lo  matched  filtering  will  also  be  treated  as  narrowband.  In  addition,  unless  stated  otherwise, 
the  analog  components  w  ill  be  assumed  to  have  a  uniform  response  over  the  system  bandwidth.  The  perti¬ 
nent  Radar  system  parameters  are  defined  in  Table  2-2. 

In  general,  the  sampled  data  obtained  from  the  column  sensors  will  be  represented  as  N^xl  spatial 
snapshot  vectors,  x.  When  space-time  processing  is  intended,  the  spatial  snapshots,  sampled  from  PRI-to- 
PRI.  are  stacked  into  (Nc  Ni)xl  space-time  snapshots,  X-  A  space-time  vector  is  created  by  stacking  (or  ras¬ 
terizing)  spatial  snapshot  vectors,  separated  in  time  by  integer  multiples  of  a  PRI.  Double-domain  matrices 
(i.e.  rasterized  vector  and  it’s  corresponding  covariance  matrix  corresponding  to  two  different  processing 
domains)  will  be  distinguished  from  single-domain  matrices  by  using  an  underline.  During  a  single  PRI, 
different  spatial  snapshots  correspond  to  different  ranges  and  therefore  will  be  referred  to  as  range-domain 
spatial  snapshots  (or  space-range  vectors  when  stacked)  in  contrast  with  the  temporal-domain  snapshots  (or 
space-time  vectors).  This  study  will  be  primarily  concerned  with  space  and  space-time  vectors. 

The  task  of  the  signal  processor  is  to  ascertain  if  a  target  with  a  prespecified  angle-of-arrival 
( AOA)  and  Doppler  frequency  is  present  in  the  data  collected  over  a  CPI  for  each  range  bin.  Thus  a  binary 
hypothesis  test  is  performed  by  the  STAP  algorithm  block  which  then  produces  detection  reports.  Under 
the  target-absent  hypothesis,  Hq,  the  data  contains  receiver  noise  and  interference  consisting  of  clutter  and 
wideband  noise  jammers; 


HQ:x  =  c-t-4  +  n  .  (2-3) 

Under  the  target-present  hypothesis,  H|,  the  data  also  contains  the  desired  target; 

H,  :x  =  as-»-c-t-j-i-n  (2-4) 

where  a  is  the  target’s  unknown  complex  amplitude.  The  specific  form  of  these  data  vectors  for  the  various 
source  types  used  in  the  numerical  analysis/simulation  results  is  the  subject  of  Section  2.1.2  thru  Section 

2. 1 .6.  A  major  goal  of  these  sections  is  to  establish  1  **  and  2"**  order  statistical  models  for  each  source  type. 
The  covariance  matrix  of  the  total  space-time  Interference  (i.e.  clutter-plus-jammers-plus-noise)  will  be 
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denoted  by  fi. 

2.1.2  Target  Model 

The  development  of  the  target  model  closely  parallels  that  given  in  [9]  and  therefore,  only  a  brief 
synopsis  is  given  here.  The  received  target  signal  is  defined  as  the  transmission  echo  from  a  moving  point 
scatterer  which  is  to  be  detected.  It  is  characterized  by  it’s  range,  angle,  and  velocity  relative  to  the  radar, 
and  the  normalized  amplitude  (i.e.  relative  to  the  ’rms’  sensor  noise)  of  the  echo  received  at  a  sensor  col¬ 
umn  output.  The  pre-envelope  of  the  transmitted  uniform  train  of  coherent  pulses  corresponding  to  a  coher¬ 
ent  processing  interval  (CPI)  is 


i{2nfo  +  <p)  ' 

s  =  a|U(t)e  .  where  u(t)=  X  “  "i^r^  (-'5) 

"•  =  0 

The  amplitude  of  a  single  pulse  waveform,  Up,  is  normalized  to  unit  energy.  Thus  the  energy  transmitted  in 
a  single  pulse  is  Ep=a,*.  The  pre-envelope  of  the  target  echo  received  by  the  n**’  sensor  is  written  as 


sfn^)  =  a^uft  -t- Tj  -  T^)e 


■<p  2v 

e  ^  ,  where  f .  =  y~ 


(2-6) 


with  received  target  amplitude  and  Doppler  a^  and  f^  respectively  (v,  is  the  radial  velocity  of  the  target  rel¬ 
ative  to  the  platform).  The  delay  to  the  Hj**’  sensor  is  composed  of  the  round  trip  delay,  tp  to  the  first  sensor 
(i.e.  the  array  reference)  and  the  relative  delay  from  the  first  sensor  to  the  nj*  sensor,  t,.  which  are  given  by 


"r"  c 


—  and 


r(4.  .0  )  -d(n  ) 


(2-7) 


=  • -sin4»cos0 


Since  the  transmitted  signal  is  assumed  narrowband,  the  relative  sensor  delay  is  presumed  insignificant  in 
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the  complex  envelope  term  of  Equation  2-6.  This  is  tantamount  to  making  the  following  approximation. 


i2nl  (i  +  i  -  I  )  i2nl.(t-T  ) 

u(t-hT^-T^)e  =u(t-T^)e 


(2-8) 


which  ignores  the  variation  of  the  received  pulse  from  sensor  to  sensor.  The  received  signal  pre-envelope 
for  the  sensor  then  becomes 


sfn^)  = 


a^e 


Kp 


u(t  -  T^)e 


i2nf^(t- 


V  '2Kfo('  +  S’ 

e 


(2-9) 


where  a  constant  phase  term  has  been  absorbed  into  the  random  phase,  (p.  It  will  be  convenient  to  use  the 
following  definitions  for  ‘normalized  spatial  frequency’,  fj,  and  ‘normalized  Doppler  frequency’,  f,; 


f 


s 


d  .  ^  ^  ,  2v,T^ 

- =  ^sin^cosB  and  f,  =  — : — 

"s  ^0  ^0 


(2-10) 


After  the  signal  has  been  processed  by  the  receiver  (i.e.  down-conversion,  matched  filtering,  and  A/D  con¬ 
version)  the  data  can  be  represented  in  complex  envelope  form.  Consider  only  those  samples  corresponding 

to  the  target  range  gate  at  limes  t^T^+nj-jp  where  n(=0 . Nj-I .  Then  the  complex  target  samples  are  given 

by 


i2jif  n  i2nf  n 

s(n^,n^)  =  a(<^.  0.  n^,  n^)e  c  ,  n^  =  0.  -  I  and  =  0 . «  1  (2-1 1) 

where  the  echo  amplitude  and  unknown  phase  have  been  combined  into  the  complex  amplitude,  a,  and  the 
following  simplifying  assumptions  have  been  made;  (I)  the  transmit  waveform  is  insensitive  to  Doppler 
shift  such  that  it  s  ambiguity  function  at  zero  delay  is  unity  for  all  Doppler  frequencies  of  interest,  and  (2) 
no  target  range  ambiguities  exist.  A  relation  between  the  Radar  system  parameters  and  the  signal-to-noise 
ratio  (SNR)  can  be  developed  (c.g.  see  [9]),  but  in  this  paper  SNR  will  be  specified  directly  instead.  It  is 
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formally  defined  as 


SNRh 


E{|a|}  PJpG,(4..e)|g((}).e)| 

a]  (4n)^NoL^R;‘ 


(2-12) 


where  is  the  system  noise  power  from  one  sensor  channel.  The  second  relation  involving  the  Radar  sys¬ 
tem  parameters,  defined  in  Table  2-2  which  assumes  identical  sensors  and  a  target  cross-section  which  is 
fixed  over  a  CPI,  is  included  for  reference.  Next,  it  will  be  useful  to  define  generic  spatial  and  temporal 
steering  vectors  as 

i2nf  i2)tf(N-l)T 

s,(fs)=[l.e  . e  ]  and 

,  (2-13) 

i2xf,  i2itf(N-l)T 

s,(f,)=[l.e  . e  I 


Note  that  fora  given  set  of  Radar  parameters,  the  complex  amplitude  of  Equation  2-1 1  is,  in  general,  a 
function  of  azimuth,  elevation  and  range.  The  case  considered  in  this  paper  assumes  each  sensor  has  an 
identical  radiation  pattern  so  that  for  each  sensor  the  complex  amplitude  is  the  same  and,  in  addition,  the 
amplitude  is  constant  over  a  CPI.  This  allows  the  target  data  samples  to  be  written  as  the  product  of  a  space- 
time  vector,  s.  and  a  complex  amplitude,  tt,  using  the  steering  vectors  of  Equation  2-1 3  as  follows: 

[s(nj.n,)]  =  a((J).  0,  R)  (s,^  ®  s  as  (f,^.  f^^)  (2-14) 

where  ‘®’  is  the  matrix  Kronecker  product,  represents  a  spatial  steering  vector  with  spatial  frequency. 
f,j,  corresponding  to  the  target’s  angular  position,  and  is  the  temporal  steering  with  normalized  target 
Doppler,  fjj. 

2.U  Receiver  Noise 

The  predominant  system  noise  component  in  a  channel  is  assumed  to  be  internally  generated  with 
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a  level  established  by  the  sensor  front-end  RF  components.  The  noise  in  each  channel  i.s  assumed  to  be 
mutually  uncorrelaied.  identically  distributed,  and  be  modeled  as  a  zero-mean  Gaussian  random  proces.s. 
For  the  case  of  the  space-time  data  vector  (i.e.  data  collected  for  a  particular  range  cell),  the  system  band¬ 
width  is  presumed  sufficiently  larger  than  the  PRF  such  that  temporal  samples  are  also  unconelated.  Under 
these  conditions  the  correlation  matrix  of  the  system  noise  is  given  by 

E„HE{nn"}  =  .  (2-15) 

For  convenience,  the  channel  noise  power  will  be  normalized  to  unity  without  loss  of  generality,  since  the 
results  in  this  paper  are  characterized  in  terms  of  noise  ratios  (e.g.  SNR,  CNR,  JNR,  etc.).  If  desired,  the 

noise  power  can  be  linked  to  the  Radar  system  parameters  by  On^=NoB.  As  previously  noted,  the  system 
noise  level  is  established  in  the  sensor  front-end,  thus  any  subsequent  linear  operations  on  the  data  (e.g. 
beamforming,  filtering,  etc.)  will  not  produce  any  significant  change  in  the  noise  ratios  (i.e.  SNR,  CNR.  and 
JNR)  specified  at  the  sensor  outputs. 

2.1,4  Clutter  Modeling 

The  main  focus  of  this  paper  is  the  study  of  techniques  to  effectively  handle  jammers  in  a  typical 
WAS  radar  interference  environment.  The  intent  of  clutter  modeling  is  to  provide  a  representative  clutter 
component  so  that  the  effects  of  jammer  processing  on  STAR  based  clutter  suppression  can  be  assessed.  For 
an  airborne  surveillance  radar  the  predominant  clutter  return  is  from  the  earth  and  this  will  be  the  only  clut¬ 
ter  component  used  for  numerical  analysis  (i.e.  simulation  results)  in  this  paper.  A  thorough  treatment  of 
clutter  modeling  typically  used  in  airborne  radar  is  provided  in  (9).  Here,  the  primary  aspects  of  the  clutter 
model  are  reviewed.  The  radar/earth  geometry  is  treated  using  the  ‘4/3  effective  radius  spherical  earth’ 
model.  Clutter  from  a  specific  range  cell  is  represented  as  the  sum  of  individual  scatters  all  at  the  same 
range  and  uniformly  spaced  In  azimuth.  The  total  clutter  component  in  a  range  gate  is  the  sum  of  contribu¬ 
tions  from  the  unambiguous  range  cell  and  any  ambiguous  range  cells.  The  signal  model  for  the  radar  echo 
from  a  clutter  patch  is  treated  in  a  fashion  similar  to  Section  2. 1 .2  for  a  target.  The  clutter  amplitudes  on 
each  sensor  from  a  single  pulse  are  identical  due  to  the  assumption  of  identical  sensor  patterns  and  the  nar- 
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rowband  signal  approximation.  However,  pulse-to-pulse  fluctuations  due  to  various  system  and  environ¬ 
mental  factors,  typically  referred  to  as  intrinsic  clutter  motion  (ICM),  will  be  modeled  as  a  wide-sen.se 
stationary  random  process  with  a  Gaussian  power  spectrum.  Using  the  spatial  and  temporal  steering  vec¬ 
tors  from  Equation  2- 1 3.  the  composite  space-time  clutter  vector  can  be  written  as 


N  N 

r  ac 


I  =  I  k  =  I 

r(4»..e,)  d 

=  — V— 


where 


2V«j)k,e,)  v 


(2-16) 


with  the  spatial  and  temporal  frequencies  of  the  (l,k)'^  clutter  patch  given  by  fj^  and  fj^-  The  large  product 
symbol,  *•*  represents  a  Hadamard  matrix  product  (i.e.  element  by  element).  The  index  T  enumerate.s  the 
ambiguous  range  bins  up  to  and  ‘k’  enumerates  the  azimuth  patches  up  to  Nj^.  The  vector  of  complex 

amplitudes,  an-,  represent  the  rentms  for  each  pulse  in  a  CPI  for  the  (Ik)'*’  clutter  patch.  The  clutter  patch 
amplitudes  are  modeled  as  wide-sense  stationary  zero-mean  random  variables  which  are  uncorrelated  from 
patch-to-patch  and  whose  variance  depends  on  various  Radar  system  parameters  (similar  to  a  target  as 
shown  in  Equation  2- 1 2  as  well  as  the  clutter  patch  position  and  reflectivity,  Ooc(4).0).  For  the  simulations 
in  this  paper  a  reflectivity  model  proposed  by  Barton  in  [  1 1  ]  for  various  types  of  land  clutter  is  used.  The 
correlation  matrix  for  the  clutter  space-time  data  vector  can  be  written  as 

Nr  N.C 

Bc-E{cc»}  =  XZ 

i  =  I 


H 

s  s 
sc  sc 


(2-17) 


As  indicated  previously,  the  temporal  fluctuations  of  the  clutter  patch  returns  are  modeled  as  having  a 
Gaussian  shaped  power  spectrum.  Hence,  the  temporal  autocorrelation  function  for  a  clutter  patch  is  given 
by 


II 


(MS) 


r^(n,)  =  exp(-2(7tai^^n,)^l 

where  the  clutter  spectral  spread  parameter,  o,cn,,  is  assumed  to  be  the  same  for  all  clutter  patches.  The  ICM 
sundard  deviation.  Ojcm.  is  the  rms  clutter  spectral  spread  normalized  by  the  system  PRF.  The  clutter  patch 
temporal  correlation  matrix  becomes 

E  n.A^  where  =  Toeplitz(  [r^.  (0) ,  ( 1 (N^  -  I )  | )  (2-19) 


and  ij,  is  the  power  in  the  (Ik)*^  clutter  patch.  The  function  Toeplitz()  represents  a  Heimitian  Toepiitz 

matrix  whose  first  row  is  given  in  the  parentheses. 

The  final  step  in  the  clutter  modeling  is  to  characterize  the  clutter  power  level.  As  previously  noted, 
the  clutter  power  level  will  be  specified  in  terms  of  CNR.  Define  the  average  total  clutter  power  per  pulse 

per  sensor*,  and  the  CNR  as 


2 


N  N  , 

N  N.  2^  2-i 

*  ‘  1  =  Ik  =  I 


11^  and  CNR  s 


(2-20) 


where  tr(]  denotes  matrix  trace.  Using  this  definition  of  clutter  power  and  a  relationship  given  in  [9)  for  the 
CNR  of  an  individual  clutter  patch,  the  space-time  correlation  matrix  can  be  written  in  a  way  which  eluci¬ 
dates  the  dependence  of  the  correlation  structure  on  the  antenna  transmit  pattern,  the  receive  sensor  pat¬ 
terns.  and  the  ground  reflectivity. 


I .  Since  an  ideal  uniform  linear  array  has  been  used  along  with  an  uncorrelated  identically  distributed  clutter  patch 
model,  the  conelation  matrix  is  doubly  block  Toepiitz  and  an  average  is  not  needed,  however  to  include  the  ca.scs 
where  the  Toepiitz  structure  is  not  preserved  this  definition  will  be  used. 
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(2-21) 


R 

“C 


r  ac 

=  z  z 


*  H, 


H 

'  s  s 
sc  sc 


where 


(G,  (4)^.  e,)  |g  (4)^.  e,)  |\,(4>k.  0,)sec  V,)  /R,^ 

Nr  N,, 

I  I  (G,(4)3.e,)|g«fr,.e,)|\,(4*3.0c)5^Ve)/Re 

a  =  I  c  =  I 


2.1^  Wideband  Noise  Jammer  (W'NJ)  Modeling 

In  this  section,  continuous  wideband  independent  noise  jammers  are  considered.  By  definition,  the 
jammer  spectrum  is  flat  and  encompasses  the  entire  receiver  bandwidth.  In  the  literature  (e.g.  (9)),  jammer 
space-time  models  for  narrowband  systems  often  assume  that  the  data  samples  are  uncorrelated  pulse-to- 
pulse  but  perfectly  correlated  sensor-to-sensor  (i.e.  across  the  array).  In  this  paper,  jammer  sample  correla¬ 
tion  will  be  modeled  based  on  a  Gaussian  shaped  receiver  spectral  characterisiic.This  will  lead  to  a  Gauss¬ 
ian  shaped  autocorrelation  function  similar  to  the  clutter  ICM  model.  The  autocorrelation  function, 
assuming  wide-sensc  staiionarity,  is  a  function  of  the  delay  between  samples,  t,  and  normalizing  to  unit 
power  is  given  by 


r.(x)  =expt^Ji(B  T)^I  .  (2-22) 

where  the  bandwidth,  B,  used  here  is  the  equivalent  noise  bandwidth.  The  pulse-to-pulse  correlation  for 
typical  PRF  and  receiver  bandwidth  combinations  (i.e.  Bt  products)  will  be  negligible,  and  therefore,  for 
space-time  data  vectors  the  jammers  will  be  treated  as  temporally  uncorrelated.  In  contrast,  the  spatial 
decorrelation  in  samples  taken  across  the  array  will  be  large  enough  to  cause  a  reduction  in  jammer  cancel¬ 
lation  and,  hence.  Equation  2-22  will  be  employed  in  the  construction  of  the  jammer  spatial  covariance 
matrix. 

Begin  by  considering  the  combined  spatial  snapshot,  j,  for  Nj  jammers  with  coordinates  (4>k.0k.Ri;) 
written  as 
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(2-23) 


Nj 

j  =  X  Jk  •  -ik  =  ^Jk  •  ^jk  • 

k  =  1 

ajjj  is  the  vector  of  random  amplitudes  received  across  the  array  of  sensors  for  (he  jammer,  and  s^j^.  is  the 
spatial  steering  vector  for  the  jammer.  For  the  case  where  the  sensors  arc  identical  (i.e,  the  case  .studied 
here),  the  elements  of  the  k*^  jammer  amplitude  vector,  are  time  delayed,  filtered,  and  sampled  versions 

of  the  k^  jammer  complex  envelope  waveform  where  (he  time  delays  are  the  sensor  delays  defined  in 
Equation  2-7.  Next,  assuming  (hat  jammer  sources  are  inde(>endent,  the  jammer  spatial  correlation  matrix 
can  be  written  as 


R,  =  E{jj"l  =e|£ 
k  =  I 
N 

r-  f  1  ^  ^  A  ^ 


k  =  t 


k  =  I 


(2-24) 


where  (he  jammer  spatial  decorrelation  is  accounted  for  by  the  matrix 


Aj^HToeplitz{  (r^(T^j^(0)),r-(T^j^(l)),  r- (N^  -  I ) )  |  } 


(2-25) 


is  the  power  received  from  the  k*  Jammer,  and  is  the  sensor-to-sensor  delay  given  by  Equation  2- 

7  for  the  k***  jammer.  The  power  level  for  each  jammer  is  specified  independently  in  terms  of  jammer-lo- 
noise  ratio  (JNR)  per  sensor  per  sample  as 


2 

JNR,. 2a  =  J,|g(vek)|'  • 


(2-26) 
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The  dependence  of  the  JNR  on  sensor  patterns  has  been  included  to  illustrate  how  they  affect  the  siaicture 
of  the  jammer  correlation  matrix.  The  constant  depends  on  jammer  p>ower  density  impinging  on  the 
array  and  several  fixed  Radar  system  parameters,  but  will  not  be  important  in  this  study.  Since  the  jammer 
samples  arc  temporally  uncorrelaled,  it  is  easy  to  see  that  the  combined  space-time  correlation  matrix  for 
the  Nj  jammers  is  given  by 


“J  N, 


N. 

L  =  I 


(2-27) 


2.1.6  Analog  Beamformer  TVansformations 

The  analog  beamformer  transformation  is  an  important  aspect  of  the  cascaded  processor  of  Figure 
2-2  and  will  have  a  significant  impact  on  the  Adaptive  Space  Processor.  It  represents  a  fixed  spatial  trans¬ 
formation  from  element  (or  sensor  channel)  space  to  a  subspace  with  equal  or  lower  dimension  which  will 
be  referred  to  as  beamspacc,  regardless  of  whether  the  transformation  develops  spatial  ‘beams*  in  the  nor¬ 
mal  sense  or  not.  In  practice  the  beamformer  acts  on  the  analog  RF  or  IF  signals,  however,  in  this  paper 
only  the  effect  on  the  digitized  baseband  signals  are  of  concern.  Unless  otherwise  stated,  the  beamformer 
will  be  assumed  to  be  adequately  represented  by  a  full-rank  matrix  operation  on  the  baseband  data  vectors 
as  indicated  in  Figure  2-3.  This  implies  that  the  beamformer’s  spectral  response  over  the  receiver  band- 


Rgure  2-3:  Analog  Beamroimer  IVansfonnation 

width  is  a  constant.  For  the  narrowband  (i.e.  <I  %  of  the  carrier)  systems  under  consideration  this  is  a  rea¬ 
sonable  assumption.  If  any  beamformers  are  frequency  dependent  they  will  be  handled  in  a  manner 
consistent  with  the  situation  at  hand.  Any  non-full-rank  beamformers  would  be  wasteful  from  a  system  cost 
perspective  and  therefore  will  not  be  studied. 
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A  specific  type  of  analog  beamfomier  will  now  be  discussed,  namely  the  Butler  Matrix,  since  it 
wijl  be  fundamental  to  the  discussion  of  the  baseline  system  architectures  as  well  as  other  architectures. 
The  Butler  matrix  is  essentially  an  analog  implementation  of  a  DFT.  It  is  a  well  known  method  for  generat¬ 
ing  a  set  of  orthogonal  beams  from  a  uniform  line  array  with  uniform  response  over  a  reasonably  wide 
bandwidth  (see  e.g.  ( 1 3]).  Consequently,  the  mathematical  description  of  the  Butler  beamformer  operating 
on  baseband  data  is  a  DFT  matrix. 

2.2  Performance  Criteria  and  Analysis  Tools 

2.2.1  Detection  Performance 

In  this  study,  detection  of  a  space-lime  target  of  known  structure  in  additive  interference  is  treated 
as  a  binary  statistical  hypothesis  testing  problem.  The  STAP  detection  algorithm  block  of  Figure  2-2  uti¬ 
lizes  a  Ney man -Pearson  approach  (see  e.g.  [12])  where  detection  probability,  Pj,  is  maximized  for  a  speci¬ 
fied  probability  of  false  alarm,  Pf.  The  binary  hypothesis  lest  under  consideration  is  expressed  by  Equation 
2-3  and  Equation  2-4.  The  target's  complex  amplitude  will  be  treated  as  a  zero-mean  complex  Gaussian 
random  variable  whose  variance  is  determined  by  the  specified  SNR  and  whose  value  is  fixed  during  a  CPI 
but  varies  CPI-io-CPI.  This  amplitude  model  is  similar  to  the  well  known  SwerlingType  I  target  model. 
Another  model  commonly  used  in  the  literature  (see  e.g.  [7],  [14],  or  [15])  treats  the  target  amplitude  as  an 
unknown  complex  constant  with  a  magnitude  commensurate  with  the  specified  SNR.  However,  the  choice 
of  target  model  will  not  significantly  affect  the  results  of  this  investigation.  Although  the  conclusions  are 
based  on  comparisons  of  different  system  architectures  and  ASP  s,  the  same  detection  algorithm  is  used  in 
each  system  and  the  impact  of  target  model  on  detection  performance  should  have  an  approximately  equal 
affect  for  each  system  studied.  The  random  amplitude  model  is  used  here  because  it  not  only  exhibits  a 
more  tractable  theoretical  detection  characteristic,  which  will  assist  in  further  developments,  but  also  is  a 
reasonable  way  to  account  for  amplitude  fluctuations  experienced  by  a  surveillance  radar  which  dwells  on 
a  target  once  per  scan. 

2.2.I-1  Optimal  Detection 

The  optimal  Neyman-Pearson  detection  statistic  is  given  by  the  likelihood  ratio  (see  e.g.  ( 1 2]).  It  is 
straightforward  to  show  that  the  optimal  detector  for  the  binary  hypothesis  test  of  Section  2.2.1  is  the  col- 
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ored-noise  maiched  filler  given  by  (see  Appendix  A) 


w 

”0 


B“*s 
s  R  $ 


(2-28) 


which  has  been  normalized  for  uniiy  variance  under  ihe  target  absent  hypothesis.  Evidently,  this  filter  turns 
out  to  be  the  same  filter  for  the  case  when  the  target  amplitude  is  a  constant  or  has  constant  magnitude  with 
uniformly  distributed  phase  (Brennan  and  Reed  [14]).  The  optimal  decision  rule  is  given  by 


H  2 

w  X 
— o  ■* 


(2-29) 


and  the  corresponding  probability  of  detection  and  false-alarm  for  the  optimal  detector  are  readily  deter¬ 
mined  as  (Park  [16]) 


=  exp  [-T]^)  and 


Pd  =  exp 


"T1  S  R  s 

•o  -  — 


/  2  H,  ^ 

(w^  [E  +  OjSS  |w^) 


-^0  1 

=  exp 

where  Oj  =  E {Ia|^}  andp^HO^s^E  *§ 


(2-30) 


The  detection  performance  is  seen  to  be  a  function  of  a  parameter  which  will  be  referred  to  as  the  general¬ 
ized  signal-io-interference-plus-noise  ratio  for  the  optimal  processor  (SINR<,  or  P^).  It  represents  an  upper 
bound  on  the  SINR  obtainable  for  the  sensor  data  being  processed. 

2.2.1.2  Generalized  Likelihood  Ratio  Algorithm  (GLR) 

The  STAR  detection  algorithm  block  of  Figure  2-2  will  utilize  a  GLR  test  for  all  the  systems  ana¬ 
lyzed.  Kelly  [  1 5]  developed  the  maximum-likelihood  GLR  algorithm  and  it’s  analytical  detection  perfor¬ 
mance  for  the  case  of  detecting  a  signal  of  known  structure  with  unknown  constant  complex  amplitude 
embedded  in  zero-mean  Gaussian  interference.  Kelly  assumed  that  a  set  of  independent  identically  distrib- 
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uted  (i.Ld.)  data  vectors  (i,e.  secondary  data)  was  available  having  the  sanie  statistics  as  target-abseni  data 
under  lest  (i.e.  primary  data  w'ithoui  a  target).  For  the  radar  problem,  the  secondary  data  is  normally 
obtained  from  range  cells  in  the  vicinity  of  but  different  from  the  primary  data,  Wang  and  Cai  [7]  extended 
the  results  to  the  case  where  the  data  is  first  processed  by  a  linear  transformation  (i.e  the  JDL-GLR).  Park 
[16]  adapted  these  results  for  the  case  of  random  complex-normal  target  amplitude  which  is  the  case  of 
interest  in  this  study.  The  results  for  a  receiver  using  Kelly’s  maximum-likelihood  GLR,  with  a  general  full- 
rank  linear  pre-transformation  applied  to  the  data,  and  employing  a  random  target  amplitude  model,  are 
summarized  below. 

The  GLR  decision  rule  for  data  which  has  been  transformed  by  a  full-rank  linear  transformation,  X 
is 


H.-l  2 

H5-I 
SjEt  • 

(  H--I  1 

d  1  -t-x-pEj  XtJ 

HI 


*oT 


(2-31) 


where  the  transform  domain  data  vectors  are  related  to  the  sensor  data  by 


-T  “  - 


X  and  Sy  =  1  s 


(2-32) 


The  unnormalized  estimate  of  the  covariance  matrix  is  the  sum  of  “K”  outer-products  of  i.i.d.  secondary 


data,  ij. 


Bt= 

k  =  I 


(2-33) 


Multivariate  normally  distributed  data  remains  multivariate  normal  under  a  full-rank  transformation  with  a 
straight-forward  alteration  of  the  mean  vector  and  covariance  matrix  (Result  4.3  in  (17)).  Thus  the  modifi¬ 
cation  of  Kelly’s  analytical  detection  results  for  transformed  data  is  straightforsv’ard.  The  transformation  I 
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used  here  represents  the  combination  of  the  linear  processing  blocks  in  Figure  2-2  preceding  the  STAR 
Detection  Algorithm  block  (i.e.  the  beamfonmer.  the  ASP,  and  the  spatial-temporal  DOF  reduction  bliKk) 
all  of  which  are  treated  as  fixed  for  the  purposes  of  analyzing  detection  performance.  The  combined  trans¬ 
formation  can  be  written  in  terms  of  the  response  matrices  of  it’s  component  blocks  as 

T=  [In^0  (LB\V)1G  (2-34) 

which  emphasizes  that  the  spatial  processing  blocks  (i.e.  L,  B.  and  W)  are  fixed  during  a  CPI.  In  this  paper, 
the  space-time  DOF  reduction  transformation,  d,  performs  temporal  reduction  only  using  a  partial  DFT 
matrix  as  discussed  in  Section  2.3. 1 .  In  this  case,  d  is  the  kronecker  product  of  selected  columns  of  a  DFT 
matrix  and  an  identity  matrix  of  size  Nj. 

From  [7],  the  probability  of  false  alarm  is  given  by 

K-N_+  I 

Pf=('-^oT)  <2-35) 

where  Nj  is  the  dimension  of  the  transformed  data  vectors,  and  the  “T’  subscript  denotes  dependence  on 

the  transformation.  From  [16],  the  probability  of  detection  for  the  random  target  amplitude  case  is  calcu¬ 
lated  as 
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f  _ K! _ 

(K-N^+ l)!(N^-2)!' 


K  -  ♦  I 


N, 


(1  -  P) 


(2-3SI 


and  the  SINR  parameter  for  transformed  data  is 


n  _  2  H„  -I 

Px-^^s^T  —  T  -T’ 


(2-39) 


The  detection  threshold  is  seen  to  be  a  function  of  the  specified  false  alarm  probability,  the  number  of  sec¬ 
ondary  data  vectors,  K,  and  the  dimension  of  the  transformed  data  vectors,  Nj.  Using  the  transformation  to 
reduce  the  data  dimension  results  in  a  smaller  threshold  for  the  same  Pf  thereby  improving  P^.  The  proba¬ 
bility  of  detection  is  a  function  of  the  specified  false-alarm  probability,  Pf,  the  transformed  SINR,  Pj.  the 
number  of  secondary  data  vectors  and  the  transformed  data  dimension.  In  order  to  more  readily  see  the 
effects  that  the  transformation  has  on  Pj  thru  Pdlp-  Equation  2-37  is  simplified  in  Appendix  B  with  the  result 


pPt+  ' 


1/(K-N_+  I) 

ppy+d/pf) 


K  -  N.p  +  I 


(2-40) 


Since  P(<1  and  for  0<p<l  and  (K-Nx+I  )>0,  it  follows  that  Pj|p  is  a  monotonically  increasing  function  of 
Px  and  therefore  so  is  P^  since  fp(p)>0  and  is  not  a  function  of  Px-  With  substantially  more  work,  it  can  also 
be  demonstrated  that  P<j|p  and  Pj  are  monotonically  decreasing  functions  of  Nxand  monotonically  increas¬ 
ing  functions  of  K.  This  is  not  proved  here  because  it  is  not  essential  to  this  study.  The  point  is  that  detec¬ 
tion  performance  can  be  improved  by  using  a  linear  transformation  to  reduce  the  data  vector  dim^ion. 
However,  the  transformed  SINR,  Px,  may  consequently  be  reduced  by  the  transformation  thereby  degrad¬ 
ing  detection.  Hence,  it  is  desired  to  design  the  transformation  so  as  maximize  SINR  for  a  given  transform- 
domain  dimension  and  to  select  the  transformed  data  dimension  which  results  in  the  maximum  improve¬ 
ment  in  Pj.  Thus,  in  the  design  of  the  signal  processor  there  are  two  main  objectives  as  far  as  detection  per¬ 
formance  is  concerned,  ( I )  maximizing  SINR  for  a  given  Nx  and  (2)  determining  the  Nx  which  maximizes 
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Pd  by  providing  the  best  balance  between  SINR  loss  and  estimation  loss. 
2.2. 1.3  Detection  Probability  Based  Performance  Metrics 


In  order  to  determine  how  different  systems  compare  in  their  detection  capability,  three  inteirelated 
performance  metrics  will  be  discussed.  The  first  metric  is  based  on  the  SINR  parameter  which  may  be  plot¬ 
ted  directly  using  a  normalized  version  of  Equation  2-39  (i.e.  normalized  with  respect  to  SNR),  or  may  be 
plotted  using  the  “SINR  Loss”  relative  to  the  optimal  SINR  (P„  given  in  Equation  2-30)  defined  as 

SINR  Loss  =  P/p.^. 

The  SINR  Loss  is  a  measure  of  the  degradation  of  the  detection  capability  (via  the  SINR  parameter)  due  to 
preprocessing  the  data  (i.e.  processing  raw  sensor  data  prior  to  generating  the  detection  test  statistic).  The 
SINR  metric  is  useful  for  comparing  systems  which  have  the  same  transformed  data  dimension,  but  it  does 
not  include  other  factors  influencing  detection,  such  as  the  loss  incurred  due  to  having  to  estimate  the  inter¬ 
ference  covariance  matrix  which  is  primarily  a  function  of  K  and  N-p. 

A  second  detection-based  performance  metric,  referred  to  as  “Detection  Loss”,  which  will  include 
all  model-based  effects  on  detection  performance  (i.e.  SINR  Loss  and  Estimation  Uss)  is  defined  as  the 
additional  SNR  required  by  the  system-under-test  in  order  to  achieve  the  same  detection  probability  as  the 

optimal  Neyman-Pearson  detector  of  Section  2.2. 1 . 1 .  Numerical  evaluation  of  Detection  Loss  proceeds  by 
determining  the  ratio 


Detection  Loss  = 


SNR| 


SNRJ 


P  .  p 

gif  ~  o,  spec 


p  ~  p 

<f.opi  d.  spec 


(2-42) 


where  the  dependence  on  the  target  structure,  interference  covariance  matrix,  specified  detection  and  false- 
alarm  probabilities,  data  transformation,  number  of  secondary  data  vectors  (i.e.  K)  and  data  vector  dimen¬ 
sion  (i.e.  N-r)  is  indicated.  The  difference  in  Detection  Loss(dB)  between  two  systems  represents  the  rela¬ 
tive  SNR(dB)  requirement  to  detect  a  target  of  a  given  structure  embedded  in  interference  with  specified 
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covariance  matrix.  It  can  be  directly  related  to  the  size  of  a  target  (i.e.  the  target  radar  cross-section)  that 
one  system  can  detect  relative  to  the  other.  It  can  also  be  indirectly  related  to  detection  range  and  system 
resources  {such  as  transmit  power,  noise  figure,  transmit  sideiobe  levels,  etc.)  and  as  such  represents  an 
impoftanl  way  to  asses  the  impact  on  practical  system  issues. 

The  third  detection-based  performance  metric,  referred  to  as  “Estimation  Loss’*,  is  based  on  a  con¬ 
venient  factorization  of  the  Detection  Loss.  Using  P^j  formulas  from  Equation  2-30.  Equation  2-36,  and 
Equation  2-40  the  SINR  required  to  achieve  a  specified  P^j  can  be  obtained.  Then  rewriting  Equation  2-42 
in  terms  of  this  required  SINR  by  making  use  of  the  SrNR  relations  in  Equation  2-30.  Equation  2-.19,  and 
Equation  2-4 1 ,  the  Detection  Loss  can  be  w'ritien  as 

Detection  Loss  =  Estimation  Loss(P^  Pp  K,  Nj)  •  SINR  Loss(R,  s,  T) 

=P.  .  (2-43) 

where  Estimation  Loss  = 

This  factorization  effectively  separates  the  dependent  variables  into  two  (almost)  disjoint  groups.  The  qual¬ 
ifier  ‘almost’  is  used  here  because  the  data  dimension,  Nj.  implicitly  influences  the  SINR  Loss  through  T 
The  Estimation  Loss  is  a  measure  of  the  degradation  of  the  detection  capability  for  a  specified  Pj  and  Pf  as 
a  consequence  of  having  to  estimate  the  interference  statistics  for  data  with  dimension  Nj  using  K  i.i.d, 
sample  vectors.  These  performance  metrics  are  power  ratios  and  are  conveniently  calculated  in  decibel 
units  resulting  in 

Detection  Loss  (dB)  =  SINR  Loss  (dB)  +  Estimation  Loss  (dB) .  (2-44) 

It  should  be  noted  that  the  required  SNR  and  SINR  of  Equation  2-42  and  Equation  2-43  are  implicit  func¬ 
tions  of  the  modeling  assumptions  and  detection  algorithms  used.  That  is  why  the  P^j  subscripts  contain 

*opt’  for  optimum  and  ‘gir’  for  the  GLR  algorithm. 

2J,J1  Jammcr-to-Noisc  Ratio  (JNR)  Improvement 

The  primary  function  of  the  ASP  within  the  cascaded  system  architectures,  which  incidentally  are 


SINRL 

opi  spec 
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ihe  focus  of  this  study,  is  to  suppress  jammer  signals,  thereby  relieving  the  STAR  of  this  task.  How  well  the 
ASP  performs  this  task  is  certainly  not  the  only  determining  factor  in  the  systems  overall  performance,  but 
oftentimes  particularly  good  or  bad  detection  performance  can  be  attributed  to  ASP  jammer  rejection  capa¬ 
bility.  Thus,  as  a  figure-of-merit  for  the  ASP,  either  the  JNR  at  the  ASP  output  or  the  ASP  jammer  cancella¬ 
tion  ratio  (JCR)  will  be  used.  The  JNR  at  the  input  of  the  ASP  is  defined  as  the  total  jammer  power  to  noise 
power  ratio. 


JNR 

in 


tr[(LB)”R.(LB)| 
tr((LB)"R^(LB)  |  ’ 


(2-4.S) 


where  the  system  spatial  noise  covariance  matrix,  R„,  is  the  identity  matrix  for  these  analyses.  The  JNR  at 

the  output  of  the  ASP  for  all  channels  combined  (i.e.  the  entire  target-clutter  subspace),  is  similarly  defined 
with  the  ASP  matrix  filter,  W,  included  as 


_tr[(LBW)"R  (LBW)J 

out  ~ '  H  •  (2-46) 

tr[(LBW)  R^(LBW)| 


Also  of  interest  is  the  JNR  in  the  one-dimensional  target  subspace  at  the  output  of  the  ASP  as  well  as  the 
JNR  m  the  one-dimensional  subspaces  comprising  the  rest  of  the  target-clutter  subspace.  These  subspaces 
can  be  thought  of  as  basis  vectors  for  the  target-clutter  subspace  and  are  defined  as  the  columns  of  the  clut¬ 
ter  Beamformer,  B^,  after  being  transformed  by  the  beamsteering  matrix,  L.  The  JNR  for  the  n'*’  clutter  sub- 
space  basis  (ihc  first  is  chosen  to  be  the  target  subspace)  is  defined  as 


_  (Lb„)”(LBW)  (LBW)”r.(LBW)  (LBW)“(Lb  ) 

-  H - H - ^  (2-47) 

(Lb^)  (LBW)  (LBVV)  R„ (LB W)  (LBW) ” (Lb  ) 


where  bp  is  the  n'**  column  of  the  clutter  Beamformer,  B^.  In  general,  the  sum  of  the  individual  JNRou,(n) 
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will  not  necessarily  equal  (he  because  (he  clutter  Beamformer  columns  are  not  constrained  to  be  an 

onhonormal  basis.  Typically.  is  used  to  investigate  the  total  jammer  suppression,  while  the  individ¬ 

ual  JNRotji(n)  are  used  to  determine  jammer  suppression  along  a  given  target-clutter  basis  vector.  The  jam¬ 
mer  suppression  along  the  target  subspace  is  usually  the  most  pertinent  since  it  provides  a  measure  of  how 
much  jammer  residue  needs  to  be  suppressed  by  the  STAR  The  jammer  cancellation  ratio  is  simply  defined 
as  the  ratio  of  oulput-to-input  JNR.  For  the  JNR  definitions  of  Equation  2-45  and  Equation  2-46  total  JCR 
is  given  by 

JCR  =  JNR  /JNR  .  (248) 

oui  in  '  ' 

For  the  individual  subspace  JNR  of  Equation  2-47  a  corresponding  input  JNR  formula  should  be  used  to 
calculate  the  JCR  for  that  subspace. 

2.2J  Extraneous  Interference  Considerations 

The  performance  metrics  considered  so  far  have  involved  the  detection  of  a  target  embedded 
within  a  particular  class  of  interference  (i.e.  complex  Gaussian  distributed  noise,  clutter  and  jammer  sig¬ 
nals)  and  the  corresponding  control  of  false-alarms  from  these  expressly  characterized  interference 
sources.  Also  of  concern  is  the  rejection  of  extraneous  signals  which  are  not  included  in  the  interference 
characterization.  These  types  of  interference  are  varied  and  difficult  to  analyze  within  a  general  framework. 
A  possible  classification  of  these  sources,  which  is  closely  linked  with  the  way  the  Gaussian  proce.ssor  of 
this  study  processes  data,  is  to  distinguish  between  those  sources  which  affect  the  primary  data  only,  those 
which  affect  the  secondary  data  only,  and  those  which  contribute  to  both.  The  first  type  affects  the  proces¬ 
sor  only  thru  it’s  contribution  to  the  primary  data,  thereby  causing  erroneous  detections  due  to  the  addi¬ 
tional  power,  or  possibly  degrading  detection  by  combining  destructively  with  the  desired  signal.  The 
second  extraneous  source  type  can  also  degrade  detection  and  cause  erroneous  detections  through  their 
influence  on  the  estimated  covariance  matrices  and  the  perturbation  of  the  assumed  probability  distribu¬ 
tions.  The  third  type  may  be  of  less  concern  if  it  is  correlated  between  primary  and  secondary  data  since  the 
matched  filter  will  tend  to  suppress  such  interference.  On  the  other  hand,  if  the  primary  and  secondary  con¬ 
tributions  are  uncorrelated  for  a  type  3  source,  then  the  effects  attributed  to  both  type  I  and  2  can  be 
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present. 


2.2J.1  Filter  Response  Pattern  Definitions 

In  this  study,  only  the  rejection  of  the  first  type  of  extraneous  interference  discussed  in  Section 
2.2.3  IS  of  interest.  A  common  way  investigating  a  processor’s  susceptibility  to  interference  sources  in  the 
primary  data,  but  not  accounted  for  in  the  design  of  the  filter,  is  to  determine  the  processor’s  response  to 
these  signals  relative  to  it’s  response  to  the  desired  signal.  For  the  surveillance  radar  system  being  consid¬ 
ered.  It  has  been  assumed  that  signal  leakage  between  range  bins  is  negligible.  However,  signals  in  the 
same  range  bin  but  with  different  AOA  and/or  Doppler  than  the  desired  signal  (i.e.  assumed  target)  can  be 
present  in  the  pnmary  data  and  need  to  be  suppressed.  The  figure-of-merit  chosen  for  comparing  the 

robustness  of  different  systems  against  such  interference  will  be  referred  to  as  the  "System  Response  Pat- 
^  and  is  defined  as 


(s,  0  s ) 
^oT?T 


where  w^.j.  s  r:J.  s.^. . 


(2-49) 


the  steenng  vectors,  s,  and  s,.  of  Equation  2- 1 3  have  been  employed,  and  the ’T’  subscript  emphasizes  the 
dependence  on  the  system’s  transformations/filters  prior  to  the  STAP  detection  algorithm  block.  Note  the 
implied  dependence  on  the  assumed  target  normalized  frequencies.  The  space-time  weight  vector.  rep¬ 
resents  the  optimum  filter  for  the  transformed  data.  It  can  also  be  related  to  the  data  filter  in  the  numerator 
of  the  GLR  test  statistic  of  Equation  2-3 1  when  the  covariance  matrix  has  been  perfectly  estimated  (or 
equivalently  when  the  number  of  i.I.d.  secondaries  increases  without  bound).  In  connection  with  this  obser¬ 
vation.  Arcan  be  viewed  as  a  normalized  ‘‘asymptotic’’  filter  magnitude  response  to  sources  different  from 
the  assumed  target.  Normalization  by  the  matched  target  response  allows  a  link  to  be  developed  between 
Ao-  and  the  extraneous-signal  SNR  handling  capability  of  a  particular  system,  as  discussed  below. 

A  more  pertinent  interpretation  can  be  assigned  to  Aj  by  referring  to  Kelly’s  work  in  [18]  which  is 
related  to  the  detection  of  mismatched  signals  using  the  GLR  algorithm.  Kelly’s  analysis  utilized  a  decont- 
position  of  the  ‘‘available  SINR’’  (i.e.  ^  of  Equation  2-39  with  Sp  replaced  by  the  steering  vector  of  the  dis- 
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placed  source)  into  two  components  orthogonal  with  respect  to  the  inverse  of  the  covariance  matrix, 

The  component  parallel  to  the  desired  steering  vector  he  referred  to  as  the  SINR  available  for  detection 
(SlNR^j.j)  and  the  other  he  called  He  showed  that  detection  performance  is  influenced  by  SENR^j^^ 

in  the  same  way  as  (Jj  for  the  matched  signal  case,  but  the  effect  of  SINRjQjj  was  substantially  more  com¬ 
plicated  since  it  resulted  in  a  modification  of  the  distribution  of  the  loss  factor  (i.e.  it’s  pdf,  which  happens 
to  be  a  Bela  distribution  for  the  matched  case,  is  now  a  function  of  the  SENRi^^  as  well  as  other  parame¬ 
ters).  The  detection  of  a  mismatched  signal  was  shown  to  be  degraded  by  two  effects  related  to  the  SINR 
components:  (I )  the  reduction  of  SINR^j,  and  (2)  the  shift  of  the  “loss  factor  pdf’  towards  greater  loss  as 
SINR|oss  is  increased.  It  can  be  shown  from  Kelly’s  results  that  the  effect  of  SINRi^^j  is  diminished  as  the 
number  of  secondaries  increase.  Kelly  demonstrated  that  a  significant  shift  in  the  “loss  factor  pdf’  corre¬ 
sponded  to  substantial  differences  in  the  SINRi^^j  (on  the  order  of  5  dB  or  more),  but  it  is  expected  for  the 
systems  considered  here  that  the  relative  S[NR|q^  levels  would  not  differ  that  much  between  systems. 
However,  because  the  SINR^^j  will  typically  exhibit  small  values  (in  the  sidelobe  regions  where  extraneous 
signals  are  of  concern),  then  the  relative  SINR^ci  levels  between  systems  can  be  substantial.  For  the  system 
comparisons  in  this  study,  it  is  presumed  that  the  reduction  in  SINR^i^j  will  be  the  primary  factor  with 
respect  to  extraneous-signal  rejection.  From  this  point  of  view,  two  systems  can  be  compared  for  their 
extraneous-signal  rejection  capability  by  determining  the  relative  extraneous-signal  SNR’s  which  would 
result  in  the  same  probability  of  detection  under  target  absent  conditions.  The  ratio  of  SINR^j  for  the  two 
systems  provides  a  good  measure  (which  becomes  belter  as  the  effects  of  SINRjqjj  become  less  of  an  issue 
for  system  comparison).  Using  Kelly’s  formulas  from  [  1 8]  the  ratio  of  detection  SINR’s  can  be  written  as 

S'NR,.,(Tj)  '  SNR^  ■ 

and,  therefore,  the  “asymptotic”  system  patterns  from  Equation  2-49  can  be  used  as  a  measure  of  the  rela¬ 
tive  extraneous-signal  SNR’s  that  can  be  tolerated  by  different  systems.  More  precisely,  for  two  systems 
differentiated  by  T|  and  T2  with  the  same  and  using  the  same  number  of  secondaries,  the  ratio  of  the 
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SNR's  corresponding  equal  extraneous-signal  detection  probability  is  determined  by  setting  the  ratio  of 
Equation  2-50  to  unity.  TTius  the  difference  in  the  System  Response  Pattern  levels  indicates  the  relative 
level  of  extraneous-signal  rejection  between  the  two  systems.  It  should  be  mentioned  that  this  measure  does 
not  include  the  aforementioned  effects  of  SINR|o„  and  also  does  not  include  the  effect  of  data  dimension. 
For  the  matched-signal  case,  a  larger  data  vector  dimension  causes  the  detection  performance  to  degrade. 
Thus,  in  comparing  two  systems  for  extraneous-signal  rejection,  it  should  be  remembered  that  if  one  sys¬ 
tem  has  a  substantially  larger  data  vector  dimension  going  into  the  GLR  processor,  it  will  incur  an  addi¬ 
tional  reduction  in  extraneous-signal  rejection  capability.  The  exact  amount  is  not  easily  accessible,  but 
using  an  estimate  for  the  effective  loss  factor  such  as  the  mean  of  the  beta  distribution,  it  can  be  shown  that 
this  extra  loss  is  less  than  2  dB  if  the  number  of  secondaries  exceeds  the  data  dimension  by  at  least  a  factor 
of  3. 

Another  response  pattern  of  interest  is  that  of  the  ASP.  Specifically,  it  is  useful  to  be  able  to  visual¬ 
ize  the  perturbation  of  the  spatial  target-clutter  subspace  response  in  the  clutter  channels  due  to  the  ASP  fil¬ 
tering,  as  well  as  the  depth,  width,  and  placement  of  the  jammer  nulls.  The  ASP  applies  a  matrix  filter  to  the 
beamformer  outputs,  and  it  is  desired  to  investigate  the  response  in  the  entire  ASP  output  subspace  (i.e.  the 
target-clutter  subspace)  as  well  as  each  output  channel  individually.  For  the  combined  subspace  of  all  ASP 
output  channels,  the  “ASP  Response  Pattern"  is  defined  as 


ASR(f^;fJ  = 


(LBW)”s 


tr((LBW)“R^(LBW)) 


(2-51) 


where  the  response  has  been  normalized  to  cancel  any  scale  factor  attributable  to  the  ASP  weight  vectors, 
w(n)  and  which  gives  the  ASR  an  interpretation  of  SNR  gain  (or  loss)  for  the  direction  corresponding  to  Sj. 
Sometimes  it  is  useful  to  plot  the  response  for  each  channel  individually,  therefore  the  “ASP  Response  Pat¬ 
tern”  for  the  n***  clutter  channel  is  defined  here  as 
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(2-52) 


ASR(n.  h 


(LBw(n))  s, 


tr[  (LBM(n))”R„(LBw(n))l 


Incidentally,  the  ASP  weight  vector  is  generally  calculated  by  assuming  enough  secondary  data  vectors  are 
available  to  perfectly  estimate  the  jammer-plus-noise  covariance  matrix.  In  this  case  then,  the  “asymptotic” 
pattern  is  the  desired  result.  Although  not  directly  linked  to  detection  performance,  the  ASP  Response  Pat¬ 
tern  can  provide  assistance  in  assessing  the  effects  of  the  ASP  without  the  added  complications  of  the 
space-time  filtering  of  clutter  and  jammer  residue. 

2.2  J.2  Characterization  of  the  System  Pattern  by  a  Scalar  Performance  Metric 

The  System  Response  Patterns  given  by  Equation  2-49  are  not  always  convenient  to  use  and  con¬ 
clusions  drawn  from  them  are  subjective  in  nature.  In  order  to  judge  a  systems’  false  detection  performance 
more  efficiently,  a  scalar  metric  is  proposed  (to  be  called  the  System  Sidelobe  Level  or  SSL)  which  is  the 
average  sidelobe  level  of  the  System  Response  Pattern.  The  average  is  performed  for  a  specific  target  (i.e. 
fixed  fjj  and  f^)  over  a  range  of  normalized  spatial  and  Doppler  frequencies  which  typically  will  not 
include  the  target  region  (i.e.  the  mainlobe).  For  uniform  non-tapered  sampling,  (in  space  and  time)  an 
orthogonal  set  of  bins  generated  from  the  original  data  all  have  equal  size  in  normalized  frequency  coordi¬ 
nates.  For  this  study,  a  source  detected  anywhere  within  the  target  bin  is  considered  a  valid  detection.  Thus, 
the  target  region  is  usually  taken  to  be  the  bin  centered  on  the  assumed  target.  In  this  paper,  the  SSL  is  cal¬ 
culated  for  I  or  2  simply  connected  rectangular  regions  in  the  space-Doppler  plane.  For  instance,  it  is 
desired  to  calculate  the  SSL  over  a  rectangular  region  which  has  an  internal  rectangular  region  excluded 
(e.g.  the  mainlobe).  The  SSL  for  a  rectangular  region  with  an  excluded  internal  rectangular  region  can  be 
calculated  by  taking  the  difference  between  the  integrated  System  Responses  for  the  two  rectangular 
regions  and  dividing  by  the  area  of  the  included  region.  The  integrated  System  Response  for  a  rectangular 
region  centered  on  the  target’s  spatial-Doppler  frequency  pair  is  defined  as 


f,,*Af,/2f„.Af/2 
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Using  this  definition*  the  SSL  for  primary  and  excluded  rectangular  region  is  calculated  from 


ssL(4r,,if,.4f,..4f„;f,^.y  = 


(2-541 


The  first  SSL  metric  uses  the  entire  space-Doppler  plane  with  the  target  mainlobe  region  excluded.  The  bin 
sizes  for  the  primary  region  and  the  excluded  target  region  between  the  halfpower  points  of  the  mainlobe 
for  uniformly  sampled  untapered  data  are 


Af,  =  1 
Af.,  =  I/N, 


Af^=  1 
Afsx  =  ‘/n; 


(2-55) 


Of  course,  the  excluded  region  defined  by  Afj,  and  Afj,  can  be  made  larger  if  the  outer  sidelobe  region  is  of 
more  interest.  A  closed  form  expression  for  the  integrated  System  Response,  Aj',  (and  therefore  the  SSL) 
can  be  obtained  by  performing  the  integrations  indicated  in  Equation  2-54.  For  the  application  considered 
in  this  paper  with  uniform  spatial  and  temporal  sampling,  a  closed  form  expression  for  SSL  is  derived  in 
Appendix  C  and  repeated  here. 


=  - ii - 

»«f5T 


(2-56) 


where  the  matrices  A,  and  Aj  defined  in  Equation  C4  are  functions  of  fy.  fjj,  Af,.  and  Afj.  A  list  of  the  SSL 
metrics  used  in  this  paper  are  given  in  Table  2-1 .  SSL!  is  the  average  Sidelobe  response  for  all  spatial  and 

Table  2-1:  Region  Definitions  for  SSL  Metrics 


SSL  Metric 

Primary 
Region  Spatial 
Dimension,  Afj 

Primary 
Region  Doppler 
Dimension,  Af^ 

Excluded 
Region  Spatial 
Dimension,  Afj^ 

Excluded 
Region  Doppler 
Dimension,  Aff, 

SSL! 

1 

I 

1/Nc 

1/Nl 
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Table  2-1:  Region  Definitions  for  SSL  Metrics 


Primary 

Primary 

Excluded 

Excluded 

SSL  Metric 

Region  Spatial 

Region  Doppler 

Region  Spatial 

Region  Doppler 

Dimension,  Afj 

Dimension,  Afj 

Dimension, 

Dimension,  Af,, 

SSL2 

1 

l/Nt 

1/Nc 

l/Nt 

SSL3 

l/Nc 

'  1 

1/Nc 

l/Nt 

SSL4 

1 

l/Nt 

3/Nc 

l/Nt 

SSLS 

1/Nc 

1 

1/Nc 

3/Nt 

SSL6 

1 

l/Nt 

.5 

l/Nt 

SSL7 

1/Nc 

1 

1/Nc 

.5 

Doppler  frequencies  excluding  the  target  mainlobe  (i.e.  a  single  space-Doppler  bin  centered  on  (f^,  fjj)). 
SSL2  is  the  spatial  sidelobe  level  for  a  doppler  width  of  one  bin  and  excluding  the  target  mainlobe. 
SSL4&6  are  the  same  as  SSL2  except  the  excluded  region  is  increased  to  determine  the  sidelobe  levels  far¬ 
ther  away  from  the  target.  SSL3  is  the  Doppler  sidelobe  level  for  a  spatial  frequency  width  of  one  bin  and 
excluding  the  target  mainlobe.  SSL5&7  are  the  same  as  SS3  except  the  excluded  region  is  increased  to 
determine  the  sidelobe  levels  farther  away  from  the  target.  The  SSL  level  are  usually  displayed  together  in 
order  as  a  vector  (or  list). 

2.2.4  Interference  Scenarios  for  Numerical  Analysis 

The  bulk  of  the  numerical  analysis  in  this  study  will  be  for  the  purpose  of  comparing  the  perfor¬ 
mance  of  various  system  architectures  and  processing  algorithms.  It  is  neither  feasible,  nor  desirable  to 
compare  the  systems  for  all  types  interference,  but  instead  to  focus  on  those  scenarios  most  pertinent  to  the 
application  being  considered.  Accordingly,  the  signal  models  detailed  in  Section  2. 1  have  been  specialized 
to  the  case  of  a  airborne  surveillance  radar  utilizing  a  uniform  linear  array  of  sensors.  The  two  types  of 
interference  considered  in  this  study  are  ground  clutter  and  WNJs.  The  spatial  arrangement  of  jammers  and 
their  power  levels  will  not  be  defined  here,  but  will  be  specified  as  needed.  This  section  is  mainly  con¬ 
cerned  with  the  specification  of  the  clutter  signal  statistics  developed  at  the  output  of  the  sensors,  and  of  the 
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Radar  system  parameters  peninent  to  the  analyses.  The  Radar  system  parameters  are  listed  in  Table  2-2  and 

Table  2-2:  Radar  System  Parameters 


Symbol 

Parameter  Definition 

Default  Value 

fo 

Radar  carrier  center  frequency 

1250  MHz 

Radar  carrier  operating  wavelength 

.24  meters 

p. 

peak  transmit  power 

- 

Tp 

transmit  pulsewidlh 

- 

B 

receiver  equivalent  noise  bandwidth 

1  MHz 

PRF 

pulse  repetition  frequency  (PRF=1/T,) 

2  KHz 

G,((J).0) 

array  transmit  power  gain 

-(uniform  taper, 
no  backlobe) 

receive  sensor  pattern  (voltage) 

Omni 

mam 

system  losses  on  receive 

- 

No 

receiver  noise  power  spectral  density 

see  Section  2.1.3 

Nc 

number  of  sensor  outputs  from  array  manifold 

32 

N. 

number  of  pulses  per  CPI 

16 

d 

sensor  spacing  in  wavelengths 

.49 

Net 

number  of  elements  in  a  sensor  column 

8 

^cl 

clement  spacing  along  the  column 

.49 

the  clutter  scenario  parameters  are  given  in  Table  2-3  along  with  default  values  for  those  parameters 
required  in  the  numerical  analyses.  The  baseline  system  described  by  these  parameters  corresponds  to  a 
narrowband,  medium  PRF  radar  using  a  medium-sized  uniform  line  array. 

The  scenarios  to  be  discussed  utilize  the  default  values  of  the  radar  and  scenario  parameters  except 
as  noted.  Scenario  A  uses  all  the  default  values.  A  substantial  amount  of  ICM  has  been  selected  which 
roughly  corresponds  to  a  clutter  velocity  standard  deviation  of  1  m/s.  Also,  this  scenario  has  the  platform 
moving  without  velocity  misalignment  (i.e.  velocity  vector  aligned  with  line  array  axis)  which  is  one  of  the 
requirements  for  a  system  using  “Displaced  Phase  Center  Antenna"  processing.  In  order  to  visualize  the 
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Power  Spectral  Density  (dB) 


Table  2-3:  Clutter  Scenario  Parameters 


Symbol 

Parameter  Definition 

Default  Value 

CNR 

total  clutier-to-noise  power  ratio  per  element 

40  dB 

Oicm 

normalized  rms  clutter  spectral  standard  deviation 

.005 

hp 

platform  altitude 

9  km 

platform  velocity  vector,  Vy,  v^],  in  m/s 

[0.,  120.,  0.] 

platform  velocity  azimuth  misalignment  angle 

0° 

R.  , 

target  range 

50  km 

e,l 

target  look  direction,  [az,  el] 

[0“.  10°) 

Nac 

number  of  clutter  scattering  cells  used  in  model 

256 

characlerislics  of  Scenario  A  clulier-plus-noise,  a  3D  mesh  plot  of  the  power  spectrum  is  provided  in  Figure 


Normalized  Spatial 
Frequency 


-0.5 


Normalized  Doppler 
Frequency 


Figure  2-4:  Power  Spectrum  for  Clutter  Scenario  A 
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2-4  The  power  sp)cctrums  in  this  section  were  generated  using  the  computed  covariance  matrices  to  create 
a  truncated  version  of  the  discrete  2D  correlation  function.  A  Blackman  window  was  applied  to  the  corre- 
lation  functions  and  a  DFT  was  performed.  Clutter  Scenario  B  is  the  same  as  A  except  that  the  transmitter 
IS  steered  50°  in  azimuth.  A  gray-scale  image  plot  of  the  clutter-plus-noise  spectrum  is  given  in  Figure  2-5 


Normalized  Doppler  Frequency  Scale  (dB) 

Figure  2-5:  Power  Spectrum  for  Clutter  Scenario  B 

In  Scenario  C.  the  steering  angle  is  broadside  as  in  Scenario  A.  but  a  25°  azimuth  velocity  misalignment  is 
introduced.  Finally,  Scenario  D  combines  the  effects  of  B  and  C  with  off-broadside  steering  of  50°  and  azi¬ 
muth  velocity  misalignment  of  25°.  Gray-scale  image  plots  of  the  clutter-plus-noise  spectrums  for  Scenar¬ 
ios  C  and  D  are  given  in  Figure  2-6  and  Figure  2-7  respectively. 

23  Baseline  Systems  Used  for  Comparative  Analyses 

This  section  is  devoted  to  the  definition  and  performance  characterization  of  the  selected  ba.seline 
ssstems  used  for  comparative  analyses.  As  already  discussed,  cascaded  ASP  /  STAP  system  architectures 
arc  of  interest  which  utilize  partially  adaptive  array  techniques  for  the  spatial  pre-suppression  of  WNJ. 
From  a  practical  viewpoint,  fixed  analog  beamforming  has  been  selected  as  the  vehicle  for  implementing 
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Normalized  Doppler  Frequency  Scale  (dB) 

Figure  2-6:  Power  Spectrum  for  Clutter  Scenario  C 
the  partially  adaptive  arrays  as  opposed  to  digitizing  every  available  spatial  channel  and  performing  some 

linear  transformation  digitally.  System  configurations  will  be  distinguished  primarily  by  the  type  of  spatial 
DOF  control  and  reduction  employed.  DOF  control  techniques  refers  to  the  type  of  fixed  analog  beamform¬ 
ing  being  used.  The  DOF  reduction  techniques  considered  here  are  of  the  adaptive  spatial  type  and  are 
designed  specifically  for  suppression  of  WNJs  by  the  ASP.  Spatial  DOF  reduction  is  indicated  in  Figure  2- 
2  after  the  ASP  function  for  the  purposes  of  maintaining  generality  in  the  block  diagram,  although  this  type 
of  processing  will  not  be  considered  here.  In  this  paper,  temporal  DOF  reduction  is  accomplished  using 
localized  DFT  bins  (as  in  the  Doppler  Domain  Localized.  DDL.  approach  [8])  after  the  ASP.  Five  .system 
configurations  are  described  here  which  will  serve  to  illustrate  various  benefits  and/or  deficiencies  of  the 
Cascaded  Space-  and  Space-Time  Adaptive  Processing  (CS-STAP)  approach,  they  are: 

•Fully-Adaptive  STAP  (FA). 

•Joint  Domain  Localized  (JDL), 

•Cascaded  Fixed  Butler-JDL(CFB-JDL). 
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Figure  2-7:  Power  Spectrum  for  Clutter  Scenario  D 
•Cascaded  Symmetric  Auxiliary-JDL  (CSA-JDL), 

•Cascaded  Selective  Butler-JDL  (CSB-JDL). 


Note  that  the  word  ‘cascaded’  distinguishes  CS-STAP  systems  from  those  not  employing  spatial  pre-sup¬ 
pression  of  jammers.  The  two  non<ascaded  systems  (FA  and  JDL)  are  used  as  representative  systems  for 
the  current  approaches  to  simultaneous  jammer-and-cluiter  suppression.  The  three  cascaded  systems  were 
selected  because  of  their  potential  to  provide  a  high-performance  practical  solution.  The  FA  system  is  a 
obvious  starting  point  for  any  adaptive  system  which  is  required  to  handle  interference  with  a  complicated 
spectrum.  The  JDL  architecture  provides  a  natural  link  between  the  FA  system  and  the  selected  cascaded 
systems  as  well  as  an  excellent  representative  approach  from  the  class  of  simultaneous  panially  adaptive 
STAR  algorithms. 

A  block  diagram  of  a  generic  Adaptive  Spatial  Processor  (ASP)  is  given  in  Figure  2-8.  In  order  to 
accentuate  the  distinction  in  functionality,  the  beamformer  has  been  split  into  two  parts,  one  intended  for 
STAP  processing.  B^,  and  the  other  for  ASP  WNJ  processing.  B,.  but  of  course,  this  is  neither  a  requirc- 
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Figure  2-8:  Generic  Adaptive  Spatial  Processor  Architecture 


ment  nor  an  implementation  restriction.  The  set  of  Mj  auxiliary  channels  (or  beams)  are  sometimes  referred 
to  as  jammer  beams,  because  (heir  primary  purpose  is  to  provide  jammer  suppression.  The  set  of  N,  chan¬ 
nels  are  sometimes  referred  to  as  clutter  beams  for  similar  reasons.  Among  the  clutter  beams,  usually  at 
least  one  of  them  is  highly  directional  and  pointing  at  the  desired  target.  Sometimes  this  beam  is  distin¬ 
guished  by  referring  to  it  as  the  target  beam. 

The  ASP  has  several  functions  as  indicated  in  Figure  2-8,  three  of  which  arc  of  interest  in  this 
paper.  These  three  are  the  ‘Secondary  Data  Extraction’,  the  ‘ASP  Weight  Calculation’  and  the  ‘ASP  Weight 
Application’  functions.  ‘ASP  Weight  Calculation’  generates  a  space-domain  matrix  filter  whose  function  is 
to  suppress  jammers  while  leaving  the  target/clutter  subspace  unperturbed  as  much  as  possible,  and  which 
reduces  the  data  dimension  from  Mj+N,  to  Nj.  ‘ASP  Weight  Application’  is  merely  a  matrix  filter  opera¬ 
tion.  An  important  feature  is  that  the  weight  matrix  is  held  fixed  during  a  CPI.  The  role  of  ‘Secondary  Data 
Extraction’  is  to  obtain  clutter-free  jammer  data  to  be  used  in  ASP  weight  calculation.  This  function  is  dis- 
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cussed  funher  in  Section  3.1 . 

The  detailed  operation  or  performance  of  the  other  two  ASP  functions  are  not  important  for  the 
puq^o.se  of  this  paper,  but  it  is  assumed  that  they  are  of  high  quality.  The  ‘DownConversion  and  A/D’  can 
be  accomplished  using  a  variety  of  techniques  for  digitization  of  a  bandpass  waveform  (e.g,  analog  syn¬ 
chronous  demodulation  and  digitization,  direct  IF  sampling  and  digital  demodulation,  etc.).  The  ‘Wave¬ 
form  Matched  Filter’  shown  in  Figure  2-8  is  applied  after  digitization  and  ASP  filtering  in  order  to  illustrate 
how  receiver  cost  can  be  reduced  by  requiring  fewer  high  quality  pulse-compression  channels,  but  this  is 
not  necessary.  The  ‘Waveform  Matched  Filter’  is  a  range-domain  pulse-compression  filter  whose  range 
sidelobes  are  assumed  to  be  negligible.  Since  the  ASP  weights  are  fixed  overa  CPI,  and  therefore  have  flat 
sp>ectral  resptonse  in  the  range-domain,  they  will  not  interfere  with  pulse-compression.  If  the  ASP  weights 
were  allowed  to  vary  in  the  range-domain  then  the  interactions  with  the  matched-filter  would  need  to  be 
addressed.  For  those  interested  in  adaptive  canceller  and  pulse-compression  interactions  see  (22)  for  exam¬ 
ple.  The  pertinent  question  here  is  the  effect  on  the  ASP  weight  calculation.  Of  importance  is  the  effect  on 
the  secondary  data  covariance  matrix  if  matched  filtering  is  performed  after  as  opposed  to  before  the  ASP 
weight  calculation.  It  turns  out  that  the  difference  in  performance  is  negligible  as  long  as  the  secondary  data 
has  been  filtered  to  a  bandwidth  no  larger  than  the  receiver  bandwidth  which  is  typically  established  by  the 
DownConversion  and  A/D’  and/or  ‘Waveform  Matched  Filter’  functions.  Reasonable  amounts  of  varia¬ 
tion  in  spectral  shape  between  the  ASP  primary  and  secondary  data  will  have  little  impact  on  the  system 
performance. 

2.3.1  Baseline  System  Architectures 

Referring  to  Figure  2-8,  the  FA  system  has  no  ASP  and  the  clutter  beamformer  transformation  is 
the  identity  matrix.  There  is  no  temporal  DOF  reduction  either.  That  is,  every  column  is  digitized  and  every 
pulse  sample  from  every  column  is  used  in  the  STAP  detection  algorithm  (ix.  the  GLRT).  The  FA  system 
has  enjoyed  much  attention  in  the  literature  (e.g.  [I4J,  [15],  and  (9J)  due  to  it’s  flexibility  and  proclaimed 
robustness.  However,  a  FA  system  for  a  medium  sized  array  will  typically  require  an  impractical  number  of 
STAP  secondary  samples,  particularly  in  a  nonhomogenous  environment.  In  reference  to  the  loss  metrics 
defined  in  Section  2.2.1,  the  FA  processor  incurs  no  SINR  Loss,  but  will  usually  have  substantial  Estima¬ 
tion  Loss  (resulting  in  substantial  Detection  Loss)  unless  many  STAP  secondary  samples  are  used. 
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There  are  several  possible  implementations  of  the  JDL  system.  The  one  chosen  here  employs  a 
Butler  bcamformer  (i.e.  an  analog  implementation  of  a  DFT)  to  generate  the  adjacent  orthogonal  spatial 
DFT  beams,  while  the  temporal  DFT  bins  are  generated  digitally  so  that  a  local  region  can  be  easily 
obtained  for  all  Doppler  bins.  This  implementation  is  used  because  of  it’s  close  relationship  with  many  of 
the  cascaded  systems  studied  in  this  paper.  Other  JDL  implementations  are  theoretically  linked,  and  should 
have  comparable  performance  under  most  circumstances.  For  an  excellent  description  of  the  JDL  concept 
refer  to  [7].  Of  course,  the  JDL  processor  will  incur  both  SINR  and  Estimation  Losses,  but  the  Estimation 
Loss  will  be  substantially  less  than  that  for  the  FA  system  when  the  number  of  appropriate  STAP  secondary 
vectors  available  is  limited. 

The  CFB-JDL  system  is  a  logical  extension  of  the  JDL  to  a  cascaded  processor.  A  fixed  set  of  Nj 
Butler  outputs  are  selected  specifically  for  STAP  clutter  processing  and  another  fixed  set  of  Mj  Butler  out¬ 
puts  are  selected  for  WNJ  suppression.  The  jammer  channels  are  used  by  the  Adaptive  Space  Processor 
(ASP),  while  the  clutter  channels  with  Jamming  removed  are  used  by  the  STAP.  For  the  baseline  processor, 
the  WNJ  beams  are  selected  to  be  adjacent  to  the  central  JDL  STAP  beams.  The  number  of  jammer  beams 
selected  is  determined  by  the  number  of  WNJ’s  the  system  is  required  to  handle. 

The  CSA-JDL  system  uses  individual  sensors  (e.g.  columns),  symmetrically  located  about  and 
close  to  the  center  of  the  array,  to  perform  WNJ  suppression.  As  indicated  by  the  JDL  appendage,  central 
Butler  beams  are  used  for  clutter  suppression  in  the  JDL  STAP.  This  beamformer  configuration  is  similar  to 
a  combination  of  two  structures  proposed  in  [1]  (overlapping  subgroups  and  symmetric  auxiliaries),  but 
with  some  important  differences;  ( I )  an  orthogonal  beamformer  utilizing  the  whole  aperture  (i.e.  a  Butler 
beamformer)  is  used  instead  of  overlapped  subarrays,  and  (2)  sensors  used  by  the  orthogonal  beamformer 
are  also  used  as  auxiliaries  for  WNJ  suppression. 

Hie  fifth  baseline  system,  CSB-JDL,  is  a  modification  of  a  jammer  suppression  technique  some¬ 
times  referred  to  as  Adaptive-Adaptive  Array  Processing  (e.g.  see  [4]).  The  basic  concept  is  to  locate  dis¬ 
crete  interferers  (i.e.  WNJ’s  in  this  paper)  and  then  form  beams  pointed  in  those  directions  in  addition  to  the 
target  direction.  The  system  studied  here  takes  the  following  approach:  (1 )  a  Butler  beamformer  supplies 
all  the  necessary  beams  to  cover  the  entire  angle  space,  (2)  the  signal  and  clutter  channels  are  preselected 
for  JDL  STAP,  and  (3)  a  direction  of  arrival  algorithm  (there  are  many  to  choose  from)  is  used  to  locate  dis- 
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Crete  tnierfcrers  using  clutter  free  data  and  the  resulting  DOA  information  is  used  to  select  the  Butler  chan¬ 
nels  closest  to  the  interferer  directions.  In  this  paper,  the  DOA  algorithms  are  not  discussed  in  detail,  and  it 
is  assumed  that  all  WNJs  are  located  accurately  enough  to  select  the  proper  beam.  A  similar  approach  was 
proposed  in  (5]  for  use  in  systems  where  low  antenna  sidelobes  was  a  primary  goal  and  only  spatially  dis¬ 
crete  interference  was  considered.  The  effects  on  clutter  suppression  were  not  investigated,  nor  was  a  prac¬ 
tical  implementation  for  a  system  combining  this  jammer  suppression  technique  with  a  clutter  suppression 
algorithm. 

2.3,2  Baseline  System  Algorithm 

A  discussion  of  the  algorithms  employed  by  the  baselir>e  systems  is  appropriate  at  this  time.  The 
detection  algorithm  employed  for  all  systems  is  the  Generalized  Likelihood  Ratio  Test  (GLR)  for  which  the 
assumptions,  test  statistic  definition,  and  theoretical  detection  performance  are  described  in  Section  2.2.1. 

The  adaptive  spatial  filter  algorithm  utilized  in  the  ASP  by  the  three  baseline  cascaded  architec¬ 
tures  (the  ASP  is  used  in  cascaded  systems  only)  is  the  Sample  Matrix  Inverse  (SMI).  The  SMI  algorithm 
(e.g.  see  [6])  has  been  chosen  because  it  converges  to  the  optimum  filter  (in  the  sense  that  SINR  is  maxi¬ 
mized)  when  given  an  adequate  amount  of  i.i.d.  secondary  data.  Indeed,  unless  stated  otherwise,  it  is 
assumed  that  enough  clutter-free  secondary  data  is  available  to  the  ASP  so  that  the  ‘tme’ jammer-plus- 
noise  covariance  matrix  can  be  used  for  analysis  purposes.  Just  how  ‘enough’ clutter-free  secondary  data  is 
obtained  is  the  subject  of  Section  3.1.  Since  the  ASP  is  a  matrix  filter  with  multiple  inputs  and  outputs, 
more  needs  to  be  said  about  the  filter  details.  The  matrix  filter  could  be  designed  as  a  whole  by  determining 
the  output  subspace  for  which  the  SINR  is  maximized  subject  to  a  set  of  predefined  constraints.  Alterna¬ 
tively.  a  filter  could  be  designed  individually  for  each  of  the  outputs  each  with  it’s  own  constraint  set.  A  ver¬ 
sion  of  the  latter  technique,  referred  to  as  the  Decoupled  Canceller,  is  used  for  the  baseline  systems.  ASP 
algorithm  alternatives  are  explored  in  Section  3.2.  In  the  baseline  ASP.  all  the  (N,-t-M,)  spatial  inputs  (i.e. 
cluner  and  jammer  channels)  are  adaptively  weighted  and  N,  separately  calculated  (Nj-i-M,)xl  weight  vec¬ 
tors  arc  applied  to  the  input  data.  Making  use  of  the  ASP  data  notation  of  Figure  2-8,  the  baseline  ASP  filter 
for  the  n'*’  clutter  channel  (i.e.  the  n*^  column  of  the  ASP  weight  matrix,  W)  is  calculated  as 
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where 


(2-57) 


=  - VB)"(LbJ  .forn  =  I . 

R^^HE{Xjx“}  .  u^sl/((Lb„)"(LB)R;;‘(LB)"(Lb^)] 

and  where  b^  is  the  n'^  column  of  the  ‘clutter-beam’  portion  of  the  beamformer  matrix  (i.e.  B^).  The  ASP 
secondary  data  is  assumed  to  be  clutter-free  jammer-plus-noise  data  with  the  same  jammer  AOA  and  JNR 
characteristics  as  in  the  primary  data.  Thus,  for  the  purpose  of  analysis  the  following  relations  are  justified, 

X.  =  (LB)”(j-i-n)  and  R  =  (LB)  ”  (R.  R J  (LB)  .  (2-58) 

j  Aj  j  n 

The  Decoupled  Canceler  operates  as  a  set  of  independent  cancelers  each  with  it’s  own  constraint  based  on 
a  distortionless  response  for  the  clutter  beam  when  calculating  the  weight  vector.  An  alternative 
interpretation  is  that  each  filter  maximizes  SINK  where  the  target's  spatial  characteristic  for  the  filter  is 

understood  to  correspond  to  the  n**'  clutter  beam.  The  intent  here  is  to  minimize  the  WNJ  power  in  each 
STAP  subspace  basis  vector  (i.e.  clutter  beam)  while  keeping  the  STAP  subspace  unperturbed. 

2.4  Baseline  Performance  Studies  and  Degree-of^Freedom  Selection 

It  is  in  this  section  that  the  system  configuration  parameters  and  performance  levels  are  established 
which  will  be  used  in  latter  sections  for  quantitative  compzuisons.  In  addition,  this  section  intends  to 
enhance  the  motivation  for  using  cascaded  systems  from  a  performance  viewpoint.  The  cascaded  systems 
can  also  have  implementation  benefits,  but  due  to  their  dependence  on  many  external  factors  (e.g.  available 
technology,  finances,  overall  platform  goals,  etc.)  these  additional  benefits  are  not  explored. 

The  remaining  step  in  the  specification  of  the  baseline  systems  is  selecting  the  number  of  spatial 
and  temporal  channels.  Of  course,  this  selection  process  is  a  function  of  the  interference  characteristics  that 
the  system  is  expected  to  handle,  as  well  as  other  practical  considerations  such  as  system  cost  (which  is  par¬ 
ticularly  associated  with  spatial  DOF).  In  this  section,  performance  in  the  clutter  scenarios  of  Section  2.2.4 
with  assumed  WNJ  environment  requirements  is  the  primary  consideration.  By  definition,  the  fully  adap¬ 
tive  system  uses  all  sensor  outputs  and  all  pulses  within  a  CPI  so  that  channel  number  selection  is  not  an 
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issue. 


In  selecting  the  number  of  temporal  channels  for  systems  employing  DOF  reduction,  the  reader  is 
reminded  that  the  primary  consideration  in  this  paper  is  the  control  and  reduction  of  spatial  DOF.  Thus  the 
number  of  temporal  channels  at  the  output  of  the  Spatial-Temporal  DOF  Reduction  function  is  chosen  large 
enough  to  adequately  handle  the  temporal  spread  of  the  clutter  without  limiting  SINR  performance,  but 
small  enough  so  that  the  estimation  loss  is  kept  to  a  reasonable  level.  All  systems  (except  the  FA  or  as  spe¬ 
cifically  noted  otherwise)  will  use  five  temporal  channels,  which  provides  a  good  compromise  for  the  cases 
studied. 

The  number  of  spatial  channels  is  selected  by  adding  the  number  needed  for  clutter  proce.ssing  to 
the  number  needed  to  handle  the  maximum  number  of  WNJs  specified  in  the  assumed  system  requirement. 
The  number  of  spatial  channels  needed  for  clutter  processing  is  chosen  as  small  as  possible  to  minimize 
estimation  loss  and  control  susceptibility  to  extraneous  signals,  but  large  enough  to  provide  'adequate' 
SINR  so  that  clutter  rejection  is  not  a  limiting  factor  in  detection  performance.  To  accomplish  these  goals 
is  not  always  easy,  as  will  be  seen.  Another  reason  to  avoid  too  many  'extra'  spatial  DOF  in  the  STAR  for 
this  study  is  that  it  is  desired  to  be  able  to  determine  if  the  ASP  is  providing  satisfactory  jammer  rejection. 
Additional  spatial  STAR  DOF  can  mask  ASP  deficiencies  by  having  the  STAR  compensate  for  inadequate 
jammer  suppression. 

Two  WNJ  system  specifications  are  considered  here.  One  system  is  required  to  handle  a  few  (i.e.  4) 
jammers.  This  system,  identified  by  the  *4-Aux*  appendage,  is  exposed  to  a  WNJ  environment  containing 
the  maximum  required.  The  other  system,  identified  by  the  M  2- Aux'  appendage,  is  required  to  handle  a 
more  considerable  number  of  WNJs  (i.e.  12),  but  it  is  exposed  to  fewer  jammers  than  required.  It  is  the  lat¬ 
ter  situation  for  which  the  cascaded  systems  are  expected  to  provide  substantial  benefits  against  sidelobe 
jammers  (from  a  performance  as  well  as  implementation  viewpoint)  since  they  can  greatly  reduce  the  com¬ 
putational  load  and  secondary  data  requirement  for  the  STAR. 

The  parameters  for  the  baseline  systems  (i.e.  number  of  channels  and  secondary  data  samples) 
have  been  selected  based  on  some  preliminary  numerical  performance  studies  with  deference  to  the  consid¬ 
erations  discussed  above.  The  number  of  spatial  and  temporal  channels  used  by  the  processing  blocks  for 
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the  baseline  systems  is  summarized  in  Table  2~4,  The  issue  of  how  many  secondaries  are  used  for  Detection 
Tabic  2-4:  Number  of  Channels  for  Baseline  Systems  Processing  Blocks 


System  /  WNJ 
Capability 

ASP  Jammer 
Channels,  Mj 

STAP  Spatial 
Channels,  Nj=Np, 

STAP  Temporal 
Channels,  Np, 

FA 

0 

Nc 

N, 

JDL/4-Aux 

0 

7 

5 

CFB-JDL/4-Aux 

4 

3 

5 

CSA-JDL/4-Aux 

4 

3 

5 

CSB-JDL/4-Aux 

adapted  as  needed  up 
to  4 

3 

1 

5 

JDL/  12-Aux 

0 

15 

5 

CFB-JDL/  12-Aux 

12 

3 

5 

eSA-JDL/  12-Aux 

12 

3 

i 

5 

CSB-JDL/  12-Aux 

adapted  as  needed  up 
to  12 

3 

5 

Loss  comparisons  is  discussed  in  Section  2.4.1  in  conjunction  with  Estimation  Loss. 

2.4.1  STAP  Estimation  Loss  and  Secondary  Data  Support 

To  determine  the  Estimation  Loss  (as  well  as  the  Detection  Loss),  the  number  of  secondary  sam¬ 
ples  (i.c.  i.i.d.  range  bins)  assumed  available  for  STAP  covariance  estimation  must  be  selected.  This  num¬ 
ber  can  vary  drastically  in  practice,  but  it  has  already  been  mentioned  that  a  (severely)  nonhomogenous 
clutter  environment  is  of  interest.  Thus  a  relatively  small  number  of  secondaries  is  implied.  In  Figure  2-9, 
the  Estimation  Loss  is  plotted  versus  the  number  of  secondaries  used  and  parameterized  with  respect  to  the 
data  vector  dimensions  corresponding  to  the  systems  considered  in  this  study.  Estimation  Loss  is  a  weak 
function  of  P^j,  and  these  curves  are  useful  over  a  wide  range  of  detection  probabilities  (at  least  for  .05  < 

<  .95).  Due  to  the  factorization  of  Detection  Loss  into  SINR  Loss  and  Estimation  Loss,  the  system  engineer 
can  study  the  effects  of  different  system  architectures  and  algorithms  on  each  loss  term  independently.  For 
the  purpose  of  illustrating  Detection  Loss  in  a  nonhomogenous  environment,  the  number  of  secondary  data 
samples  used  for  STAP  covariance  estimation  is  selected  to  be  120  (unless  indicated  otherwise).  Note  that 
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this  number  is  insufficient  to  generate  an  invertible  covariance  matrix  for  the  FA  systems  considered  here. 
Therefore,  Detection  Loss  of  the  FA  systems  will  be  based  on  using  a  secondary  data  sample  set  10%  larger 
than  the  space-time  FA  data  vector  dimension  (i.e.  K=  I .  I  *Nc*N,).  The  FA  results  are  meant  for  reference 
only  ;  as  an  indication  of  the  FA  system  capability  if  the  number  of  'adequately  i.i.d.’  range  bins  available  is 
minimally  sufficient  to  perform  FA  processing. 

2.4.2  Performance  for  Systems  Required  to  Handle  a  Few  WNJs 

The  ‘4-Auxiliary’  system  is  required  to  handle  up  to  4  WNJs  for  each  of  the  clutter  scenarios,  and 
therefore  the  number  of  auxiliary  channels  is  Mj=4.  Through  numerical  analyses,  it  has  surfaced  that  for 
most  of  the  scenarios  investigated,  3  spatial  STAP  DOF  are  sufficient  for  clutter  suppression,  and  therefore 
Nj=3  spatial  DOF  is  chosen  for  the  Baseline  systems.  However,  for  the  case  with  2  clutter-region  jammers. 
5  spatial  STAP  DOF  are  needed  to  obtain  good  performance  and  hence  this  number  is  sometimes  used. 
Some  results  for  the  ‘4-Auxiliary’  case  are  given  in  the  following  subsections. 
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2.4.2. 1  Performance  in  Ihe  Absence  of  Jammers 

The  first  step  in  the  analysis  of  the  cascaded  configurations  is  to  determine  if  they  will  degrade  sys¬ 
tem  performance  when  presented  with  a  scenario  for  which  they  are  not  needed  (i.e.  WNJ  free  interfer¬ 
ence).  Using  the  clutter-only  scenarios  of  Section  2.2.4  the  five  baseline  systems  are  compared. 
Representative  Detection  Loss  results  are  given  in  Figure  2-10  and  Figure  2-11,  The  cascaded  systems  per- 


12 


10 


cn 

i/i 

t/i 

o 


o 

a 

O) 

O 


-  FA  :  JNnftfvK<«)“t  *53  I 

H  JOL(7.5)  :  JNR0nJ(<ie)-t-53I  .  , 

4-  CFB{4).J0U3.5>  :  -53  -6t  I 

O  CSA(4>-X>U3.5)  :  JNfnr»n>jtl(dB)*(  -60  •€!  | 

X  CSe<0)  JOU3.5)  ;  JNRr»'\->/tK<3e)«{  -61  -611 

FA  :  SSL(<S8)H  -3  4  9  -2  5  -6  3  | 

JOL  :  SSUdBH -2  4  1C  -2  6  -24  -8) 
CFB  SSU<Je)-(-2  5  1C  -13  5  -<2  -6 1 

eSA  ;  SSL(<58H -2  5  1C  -13  5  -42-6  1 

CSB  :  SSU<58H  *2  5  1C  -13  5-42-6  1 

Jammer  Fraq.  «  [  ] 

CNR  O  e4am  «  40  dB 

|Pta.  •  Covar  Samples]  •  '  te-05.  120  | 


01— 

-0.5 


-0  4 


-0.3 


-0.2 


0.3 


-0.1  0  0.1  0.2 
Normalized  Target  Doppler 
Figure  2-10:  Detection  Loss  without  Jammers:  4-Aux/Scenario  A 
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form  at  least  as  well  as  reference  systems  (i.e.  the  FA  and  JDL  systems)  and  exhibit  no  detrimental  effects 
for  clutier-only  scenarios.  Note  that  all  the  reduced  DOF  systems  exhibit  significantly  lower  spatial  SSL 
levels  (response  plots  not  shown)  away  from  the  mainlobe  as  indicated  by  SSL4&6  shown  in  the  figure  leg¬ 
ends.  Also  note  that  the  Detection  Loss  for  the  FA  system  is  not  a  function  of  target  Doppler,  because  it 
incurs  no  SINR  Loss  and  Estimation  Loss  is  independent  of  5,  R,  and  T. 

2.4^.2  Performance  Studies  for  Sidelobc  WNJ  Environments 

WNJs  located  in  the  “sidelobc  region”  facilitate  the  analysis  of  jammer  suppression  without  the 
added  complication  of  clutter  suppression  degradation.  Stated  another  way.  sidelobc  jammer  suppression 
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Figure  2-11:  Detection  Loss  without  Jammers:  4-Aux/Scenario  D 


will  have  much  less  effect  on  STAP  performance  than  clutter  region  jammers  as  long  as  the  number  of  jam¬ 


mers  does  not  exceed  the  number  of  auxiliary  channels.  The  clutter  scenarios  used  in  Section  2.2.4  are  used 


along  with  sidelobe  WNJs  to  compare  the  five  baseline  systems.  The  JNRs  for  each  jammer  will  always  be 


the  same  (unless  stated  otherwise)  and  for  these  results  are  set  to  60dB.  Detection  Loss  results  are  given  in 


Figure  2- 1 2  and  Figure  2- 1 3  for  clutter  scenarios  A  and  D  respectively.  System  Response  spatial  patterns 
are  given  in  Figure  2-14  and  Figure  2-15  for  a  system  nined  to  a  target  with  the  spatial  and  Doppler  fre¬ 
quencies  specified  in  the  plot  legend.  The  System  Response  spatial  pattern  plot  is  a  spatial  cut  from  the 
two-dimensional  System  Response  at  the  specified  target  Doppler.  Note  that  theCSA  system  produces  SSL 
levels  in-between  that  of  the  FA  system  and  the  other  systems  and  also  exhibits  more  detection  loss  than  the 
other  cascaded  systems  for  target  Dopplers  approaching  the  mainlobe  clutter.  So  far  Detection  Loss  has 
been  used  for  performance  comparisons.  From  this  point  on  it  is  desirable  to  focus  in  on  SINR  Loss  since 
the  system  architecture  and  algorithm  differences  will  be  primarily  a  function  of  5.  K.  and  I.  If  Detection 
Loss  is  of  interest  for  a  particular  level  of  secondary  sample  support  then  the  appropriate  Estimation  Loss 
from  Figure  2-9  can  be  added  to  the  SINR  Loss  results.  Thus,  in  Figure  2-16  and  Figure  2-17  SINR  Loss 
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Figure  2-12:  Detection  Loss  with  Sidelobe  Region  Jammers:  4-Aux/Scenario  A 
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Figure  2-13:  Detection  Loss  with  Sidelobe  Region  Jammers:  4-Aux/Scenario  D 
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Figure  2-14:  Spatial  System  Response  with  Sidelobe  Region  Jammers:  4-Aux/Scenano  A 
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results  are  included  for  completeness  and  ease  of  comparison  with  future  results. 

1.422  Performance  Studies  with  WNJ’s  in  the  Clutter  Region 

As  mentioned  previously,  it  is  presumed  that  the  clutter  subspace  defined  by  the  Nj  spatial  STAP 
channels  has  been  designed  in  order  to  provide  good  clutter  suppression.  It  is  desirable  to  minimize  pertur¬ 
bation  of  this  subspace  while  rejecting  WNJ  power.  However,  when  a  jammer  angle  is  located  within  the 
clutter  region  noticeable  perturbation  is  inevitable.  Hence,  it  is  of  interest  to  determine  which  system  archi¬ 
tectures  are  robust  under  these  circumstances.  Scenario  A  results  are  given  in  Figure  2-18  and  Figure  2-19 


Figure  2-18:  SINR  Loss  with  Gutter  Region  Jammers:  4-Aux/Sccnario  A 


with  two  of  the  four  jammers  located  within  the  clutter  region,  but  not  in  the  mainlobe  (i.e.  target  region). 
Significant  SINR  loss  is  evident  for  the  CSA-JDL  system  possibly  indicating  a  fundamental  problem  with 
this  approach  for  clutter  region  jamnoers.  The  other  two  cascaded  systems  showing  performartce  compara¬ 
ble  to  the  JDL.  Reduction  in  SINR  for  all  systems  (including  the  maximum  attainable  given  by  the  FA  sys¬ 
tem)  is  evident  by  comparing  the  System  Response  of  Figure  2-19  to  Figure  2-14  at  fj=0.  Note  that  all  the 
reduced  DOF  systems  exhibit  more  SINR  Loss  than  in  the  sidelobe  jammer  case.  Also  note  that  the  output 
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jammer  levels  (shown  in  temis  of  JNR  in  the  plot  legends)  are  not  suppressed  to  the  noise  floor  as  in  previ¬ 
ous  results.  Scenario  D  performance  results  with  two  clutter  region  Jammers  are  given  in  Figure  2-20.  The 
results  are  similar  to  those  for  Scenario  A  except  that  the  CSB-JDL  system  is  showing  slightly  poorer  per¬ 
formance.  Substantial  performance  improvement  can  be  achieved  by  increasing  the  number  of  spatial 
STAP  channels  to  N5=5  as  shown  in  Figure  2-21  and  Figure  2-22.  Evidently,  having  two  jammers  in  the 
clutter  region  requires  more  STAP  DOF  in  order  to  maintain  a  low  SINR  loss.  Increasing  the  number  of 
ASP  DOF  (i.e.  Mj).  however,  does  not  provide  the  same  performance  improvement. 

2.4,2.4  Performance  Studies  with  WNJ’s  in  the  Target  Region 

The  most  severe  angular  location  for  a  jammer  is  near  the  desired  signal.  In  this  case  it  is  of  interest 
to  determine  if  the  cascaded  approach  to  WNJ  cancellation  is  still  viable,  or  if  the  jammer  now  needs  to  be 
handled  at  the  same  time  as  the  clutter.  Said  another  way,  when  the  jammer  resides  within  the  target  .sub- 
space  it  may  be  advisable  to  supply  the  STAP  filter  with  all  the  interference  information.  Performance 
results  for  the  case  with  one  jammer  in  the  target  region  and  three  sidelobe  jammers  is  shown  in  Figure  2- 
23  for  Scenario  A  and  Figure  2-24  for  Scenario  D.  The  SINR  Loss  for  the  JDL  and  CFB-JDL  arc  very 
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Figure  2-20;  SINR  Loss  with  Clutter  Region  Jammers:  4-Aux/Scenario  D 


Figure  2-21:  SINR  Loss  with  Clutter  Region  jammers:  4-Aux/Nj=5/Scenario  A 
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Figure  2-22:  SINR  Loss  with  Gutter  Region  Jammers:  4-Aux/Ns=5/Scenario  D 
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Figure  2-23:  SINR  Loss  with  Target  Region  Jammer:  4-Aux/Scenario  A 


Figure  2-24:  SINR  Loss  with  Target  Region  Jammer:  4-AuxyScenario  D 
close,  while  the  CSB-JDL  system  exhibits  quite  robust  behavior  for  this  case.  A  spatial  cut  of  the  System 


Response  for  Scenario  A  is  given  in  Figure  2-25.  As  expected,  the  SINR’s  at  fj=0  are  greatly  reduced  for  all 


the  systems  including  the  SINR  Bound  (or  FA  system). 

2.4  J  Performance  for  Systems  Required  to  Handle  Many  WNJs 


The  ‘12- Auxiliary’  system  is  required  to  handle  up  to  12  WNJs  for  each  of  the  clutter  scenarios. 
Accordingly,  the  number  of  auxiliary  channels  is  set  to  Mj=12.  The  WNJ  scenarios  will  contain  four  jam¬ 
mers  as  in  Section  2.4.2,  but  this  time  the  number  of  auxiliary  channels  exceeds  the  number  of  jammers  by 
a  factor  of  three.  The  number  of  spatial  STAP  channels  is  selected  as  Nj=3  as  in  Section  2.4.2.  Selected 
results  for  the  '12- Auxiliary’  case  are  given  in  the  following  subsections. 

L4  J.l  Performance  Studies  for  Sidelobe  WNJ  Environments 


The  cascaded  systems  show  better  detection  performance  than  the  JDL  and  FA  systems  for  the 
sidelobe  jammer  case  shown  in  Figure  2-26  for  scenario  A  and  Figure  2-27  for  Scenario  D.  They  also 
exhibit  better  SSL  levels  as  indicated  in  the  plot  legends  and  as  demonstrated  in  the  System  Response  spa- 
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Normalized  Spatial  Frecjuency 

Figure  2-25:  Spatial  System  Response  with  Target  Region  Jammer:  4-Aux/Scenario  A 
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Figure  2-26:  Detection  Loss  with  Sidelobe  Jammers:  12-Aux/Scenario  A 
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Normalized  Target  Doppler 

Figure  2-27:  Detection  Loss  with  Sidelobe  Jammers:  12-Aux/Scenario  D 


(iai  cut  of  Figure  2-28  for  Scenario  A.  Similar  results  are  obtained  for  Scenario  D.  Thus,  the  cascaded  sys¬ 


tems  have  better  detection  performance  as  well  as  more  protection  against  extraneous  signals  in  the  test 


range  cell.  Since  SINR  Loss  is  used  for  comparative  analyses,  it  is  shown  here  in  Figure  2-29  and  Figure  2- 


30.  Comparing  the  SINR  Loss  with  the  Detection  Loss  reveals  a  substantial  Estimation  Loss  for  the  JDL 


system  of  about  5  dB  for  the  case  where  1 20  secondary  samples  are  used. 


2.43.2  Performance  Studies  with  WNJ’s  in  the  Clutter  Region 


As  in  Section  2.4.2.3  with  clutter  region  jammers,  it  is  expected  that  the  cascaded  systems  will 
show  some  performance  degradation,  but  now  the  JDL  should  not  degrade  as  much  since  it  has  many  extra 
DOF  to  use.  These  expectations  are  realized  as  shown  in  Figure  2-3 1  for  Scenario  A  and  in  Figure  2-32  for 
Scenario  D.  Note  that  the  CFB-JDL  system  in  particular  maintains  relatively  good  SINR  performance. 
Although  the  performance  degradation  for  WNJs  in  the  clutter  is  not  as  pronounced  as  in  the  ‘4-Aux*case. 
the  performance  of  all  the  cascaded  systems  can  be  still  be  noticeably  improved  if  the  number  of  STAR  spa¬ 
tial  channels  is  allowed  to  increase  to  Nj=5  as  evidenced  by  Figure  2-33  for  Scenario  A. 
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System  Filler  Response  (dB) 


Normalized  Spatial  Frequency 

Figure  2-28:  Spatial  System  Response  with  Sidelobe  Jammers:  12-Aux/Scenario  A 


Normalized  Target  Doppler 

Figure  2-29:  SINR  Loss  with  Sidelobe  Jammers:  12-Aux/Scenario  A 
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Figure  2-30:  SINK  Loss  with  Sidelobe  Jammers:  12-Aux/Scenario  D 
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Figure  2-32:  SINR  Loss  with  Clutter  Region  Jammers:  12-Aux/Scenario  D 
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Figure  2-33:  SINR  Loss  with  Clutter  Region  Jammers:  12-Aux/Nj=5/Scenario  A 


2.4.3.3  Performance  Studies  with  WNJ’s  in  the  Target  Region 

The  performance  of  the  CFB-JDL  and  CSB-JDL  systems  remain  competitive  with  the  JDL  for  the 
case  with  one  jammer  in  the  target  region  and  three  sidelobe  jammers  for  both  Scenario  A  and  D  as  evi¬ 
denced  by  the  SINR  Loss  results  in  Figure  2-34  and  Figure  2-35.  However,  the  CSA-JDL  system  shows 


considerable  SINR  degradation  as  was  the  case  for  the  M-Aux*  systems. 

2.4.4  Discussion  of  Baseline  Performance  Results 

The  preceding  performance  studies  have  shown  that  cascaded  systems  are  a  viable  alternative  to 
the  combined  clutter  and  WNJ  suppression  problem,  and  are  the  preferred  approach  in  circumstances  with 
limited  secondary  support  in  terms  of  Detection  Loss  as  well  as  extraneous  signal  protection.  Even  under 
circumstances  with  jammers  within  the  cluuer  and/or  signal  subspace,  for  which  the  cascaded  approach  is 
at  a  disadvantage  (since  clutter  information  is  unavailable  to  the  ASP),  the  cascaded  systems  exhibit  com¬ 
petitive  results.  In  addition,  very  little  was  said  of  the  implementation  benefits.  It  is  easy  to  see  that  reduc¬ 
ing  the  number  of  channels  that  must  be  well  matched  and  digitized  can  greatly  improve  system 
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Figure  2-35:  SINR  Loss  with  Target  Region  Jammers:  12-Aux/Scenario  D 
affordability  and  reliability.  The  cascaded  systems  can  also;  ( 1 )  reduce  the  number  of  waveform  matched 

filters  (i.e.  pulse  compressors)  which  lessens  the  computational  load,  (2)  reduce  the  number  of  STAR  chan¬ 
nels  thereby  preserving  the  performance  of  adaptive  systems  in  severely  nonhomogenous  interference 
environments,  and  (3)  diminish  System  Response  sidelobe  levels  by  eliminating  superfluous  DOF  in  the 

STAP. 
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3  Implementations  and  Alternatives  for  the  Adaptive  Space  Processor 

In  ihis  section,  attention  is  focused  on  the  implementation  of  the  ASP  baseline  approaches.  The  pri¬ 
mary  topic  is  how  to  obtain  adequate  clutter-free  secondary  data  for  WNJ  suppression.  In  the  past,  clutter- 
free  secondary  data  extraction  has  been  performed  by  using  data  which  is  known  to  have  none  or  little 
amounts  of  clutter  power.  Three  mainstream  techniques  are  mentioned  here.  Low-clutter  range  cells  can 
sometimes  be  collected  during  a  Low  PRF  mode,  but  are  restricted  to  only  a  portion  of  a  PRI.  Low-clutter 
Doppler  regions  can  be  obtained  during  a  high  PRF  mode,  but  require  special  Doppler  filtering  to  isolate 
low-clutter  regions.  This  technique  is  usable  primarily  for  circumstances  resulting  in  an  uncomplicated 
space-time  clutter  spectrum.  Clutter-free  collection  during  a  passive  mode  (i.e.  non-transmitttng)  of  the 
Radar  is  possible,  but  severely  restricts  the  time  period  for  collection.  These  techniques  all  have  restrictions 
on  when  and  how  much  clutter-free  data  is  available.  It  is  desired  to  have  little  or  no  restrictions  on  WNJ 
data  collection  so  that  an  entire  RPI  data  set  can  be  used  if  necessary.  In  this  way.  Estimation  Loss  for  WNJ 
covariance  estimation  can  be  made  negligible.  Additional  benefits  to  unrestricted  secondary  collection  are 
(1)  robustness  against  blinking  jammers,  (2)  very  fast  adaptation  possible  (even  sub-PRI  in  some  cases), 
and  (3)  provides  a  system  foundation  amenable  to  handling  the  jammer  multipath  (i.e.  hot  clutter)  problem. 
Section  3.1  presents  an  approach  which  can  realize  these  advantages. 

3.1  Sideband  Canceller  (SBC)  Concept 

This  section  develops  a  new  technique  for  generating  the  adaptive  spatial  weights  used  to  suppress 
WNJ  for  a  general  CS-STAP  system  architecture.  The  central  idea  is  to  obtain  clutter  free  WNJ  data  using 
a  narrow  RF  sideband  which  is  adjacent  to.  but  disjoint  from,  the  Radar  receiver  mainband.  The  assumption 
is  that  the  spectrum  of  a  WNJ  contributing  to  the  Radar  receiver  signal  will  also  encompass  a  non-zero 
bandwidth  surrounding  the  receiver  band.  An  additional  assumption  is  that  an  IF  or  baseband  filter  can  pro¬ 
vide  a  narrow  WNJ  passband  adjacent  to  the  Radar  band,  but  free  from  any  Radar  effected  signals  (i.e.  clut¬ 
ter  and  targets).  A  conceptual  block  diagram  of  the  SBC  is  shown  in  Figure  3-1.  Note  that  the  beam.steering 
and  beamforming  transformations  are  now  considered  to  be  functions  of  the  mainband  and  sideband  center 
frequencies.  Thus,  where  applicable,  the  response  of  the  sideband  channels  for  different  carrier  frequencies 
is  accounted  for  in  the  analyses.  The  design  of  the  sideband  canceler  involves  the  definition  of  the  Sideband 
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Figure  3-1:  Sideband  Canceler  Conceptual  Block  Diagram 
Beamsteering.  L,(f,).  and  Sideband  Beamformer.  which  may  or  may  not  be  the  same  type  as  the 

mainband  beamformer  B(fo).  the  specification  of  the  Sideband  Filter,  h,.  and  the  design  of  the  ASP  Weight 
Calculation  Algorithm. 

The  beamsteering  transformation  response.  L(f).  for  the  different  earner  frequencies  can  be  a  fac¬ 
tor  in  the  performance  of  the  SBC  concept,  particularly  when  the  Radar  look  direction  is  other  than  broad¬ 
side.  Since  a  phased-array  is  considered  here,  beamsteering  takes  the  form  of  progressive  sensor-to-.sensor 
phase  shifts.  Array  calibration  can  be  included  in  the  beamsteering  transformation  model,  but  will  not  be 
included  here.  Hence.  L(0  is  taken  to  be  a  unitary  diagonal  matrix  whose  diagonal  elements  may  be  func¬ 
tions  of  the  carrier  frequencies.  However,  sensor  response  variations  over  the  bandwidth  of  the  receiver’s 
mainband  (or  sideband)  will  be  assumed  negligible  for  all  system  components  as  discussed  in  the  System 
Model  descriptions  of  Section  2. 1 . 

Several  SBC  approaches  will  be  considered  in  this  section.  Although  it  is  important  to  keep  in  mind 
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ihc  dependence  of  performance,  processing  complexity,  and  cost  on  the  type  of  Sideband  Beamformer 
selected,  performance  is  the  primary  focus  of  this  investigation.  The  specific  Sideband  Beamformer  defini¬ 
tion  for  each  approach  is  discussed  as  needed.  Note  that  the  number  of  beams  at  the  output  of  the  Sideband 
Beamformer  is  not  necessarily  the  same  as  for  the  mainband,  Nj,.  Some  sideband  compensation  techniques 
will  require  more  beams  in  order  to  obtain  satisfactory  performance,  and  in  general  for  this  study 
With  regard  to  system  cost  it  is  desirable  to  minimize  the  number  of  receiver  channels.  However,  as  will  be 
pointed  out  later,  the  bandwidth  of  a  sideband  channel  is  kept  quite  narrow,  and  therefore  the  cost  can  be 
substantially  less  than  that  of  a  mainband  channel.  Therefore  a  large  number  of  sideband  channels  does  not 
necessarily  preclude  the  use  of  such  an  approach.  The  ASP  weight  matrix  will  always  have  columns 
(and  Nj  rows)  in  order  to  be  compatible  with  the  mainband  data  for  which  it  is  intended.  Thus,  for  sideband 
systems  with  Nf>Nb  the  ASP  Weight  Calculation  block  will,  as  part  of  the  algorithm,  provide  the  necessary 
dimension  reduction. 

The  sideband  filter  is  assumed  to  be  the  principal  factor  in  determining  the  frequency  response  for 
the  sideband  channels,  as  is  the  mainband  filter  for  the  mainband  channels.  As  discussed  in  Section  2,1 .4 
and  Section  2. 1 .5,  the  frequency  response  of  the  channels  from  the  sensors  to  the  input  of  the  STAP  are 
modeled  as  having  a  Gaussian  shaped  bandpass  for  the  purposes  of  calculating  the  data  covanance  matri¬ 
ces.  Thus,  the  exact  specification  of  the  sideband  and  mainband  filters  is  not  necessary  for  the  analyses  of 
this  study;  instead,  the  parameters  which  need  to  be  specified  are  the  center  (or  equivalently  carrier)  fre¬ 
quencies,  (q  and  f|,  and  the  corresponding  bandwidihs.  The  sideband  frequency  offset,  f^f^,  is  chosen  as 

small  as  possible  while  allowing  the  sideband  filler  to  sufficiently  suppress  radar  generated  signals.  The 
bandwidihs  arc  used  in  characterizing  the  sensor-to-sensor  decorrelation  of  the  WNJ  signals  as  discussed  in 
Section  2.1.5.  The  mainband  bandwidth  specified  in  Table  2-2  is  used  for  the  majority  of  the  results  with 
any  exceptions  specifically  noted.  The  sideband  bandwidth  is  chosen  considerably  smaller  than  the  main- 
band  bandwidth  for  two  reasons;  (1 )  minimize  the  cost  of  the  sideband  channels,  and  (2)  restrict  the  ASP  to 
respond  to  the  narrowband  spatial  characteristics  of  the  Jammers  without  consuming  spatial  DOF  to 
attempt  to  handle  bandwidth  effects.  Bandwidth  effects  are  belter  handled  using  joint  space-time  process¬ 
ing  and  therefore  arc  left  for  the  STAP. 
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The  baseline  ASP  weight  calculation  algorithm  of  Equation  2-57  needs  to  be  adapted  for  use  in  the 
SBC  concept,  in  order  to  account  for  the  characteristics  of  the  sideband  data.  Further,  the  particular  adapta¬ 
tion  technique  used  will  have  ramifications  on  performance,  cost,  and  the  type  of  sideband  beamformer 
(and  beamsteering)  used.  This  topic  is  addressed  mainly  in  Section  3.3.  but  is  also  considered  as  necessary 

in  other  sections. 

3.2  Sideband  Data  Modeling 

A  complete  description  of  the  target,  clutter  and  WNJ  data  models  for  the  Radar  mainband  was 
given  in  Section  2.1 .  Here,  the  appropriate  model  for  the  sideband  data  is  developed  in  order  to  form  a 
mathematical  basis  for  further  developments  and  analyses.  In  Section  2.1.  space-lime  data  was  arranged  by 
stacking  spatial  snapshots  for  a  given  range  bin.  producing  an  angle-Doppler  transform  domain.  In  this  sec¬ 
tion.  Doppler  processing  is  not  utilized.  Instead,  the  data  is  described  in  terms  of  spatial  snapshots  as  fime- 
tions  of  time  (discrete,  m.  or  continuous,  t)  without  any  special  regard  to  PRI  boundaries.  CPI  boundanes. 
however,  are  still  important,  because  the  ASP  weight  matrix  is  allowed  to  change  at  CPI  boundaries  only 
(i.e.  not  during  a  CPI)  for  a  particular  range  bin.  However,  the  ASP  filter  may  be  different  for  different 
range  bins.  This  constraint  on  the  ASP  weights  is  convenient  for  analysis  and  should  be  useful  for  many 
practical  situations,  but  should  not  be  considered  as  a  general  restriction  on  the  SBC  approach. 

The  Sideband  filter  is  assumed  to  have  sufficiently  narrow  bandwidth  that  the  narrowband  approx¬ 
imations  used  in  Section  2.1  are  applicable,  and  the  data  at  the  output  of  the  Sideband  filter  can  be 
expressed  in  terms  of  complex  envelopes.  Using  Equation  2-16  and  Equation  2-23  with  appropriate  modi¬ 
fications.  the  spatial  snapshot  data  model  as  a  function  of  continuous  time  at  the  output  of  the  Sideband  fil¬ 
ter  is  given  as 


x,(t)  =  (L,B,) 


(a(t)©h,(i))s^(f„(f,))+ £  £  (c,q(t)©h,(t))s^(f,^(f,)) 

l  =  lq  =  l  (3.1) 


N- 


+  X  (jk(‘)©h,(0)  +n(l)©h,(t) 

k  =  I 


where  ©  denotes  convolution,  the  dependence  on  carrier  frequency  is  indicated,  and  notation  has  been 
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slightly  modified  to  incorpxjrate  time  functions.  As  is  customary,  convolution  of  a  scalar  function  with  a 
vector  implies  convolution  with  each  element  individually.  As  mentioned  earlier,  it  is  assumed  that  the 
sideband  is  separated  far  enough  from  the  mainband  that  the  Sideband  filter  can  reduce  Radar  generated 
signals  to  negligible  levels.  Thus  after  sampling  (at  greater  than  the  Nyquist  rate),  the  discrete  sideband 
data  can  be  modeled  as 


x,(m)  =  (L,B,) 


.k=  I 


(3-2) 


where  n(m)  is  the  discrete  filtered  noise  sequence  and  ajv.(m)  is  the  discrete  filtered  jammer  sequence  for 

the  k***  jammer.  The  reader  is  reminded  that  the  beamsteering  and  beamformer  transformations.  L  and  B, 
are  in  general,  functions  of  the  carrier  frequency.  The  sideband  data  is  assumed  to  be  wide-sense  stationary 
and  a  spatial  snapshot  at  any  given  sample  has  zero-mean  and  covariance  matrix  given  by 


where  the  jammer  sequence  correlation  matrix.  A^;..  has  the  same  form  as  that  given  by  Equation  2-25, 
except  that  the  correlation  function  of  Equation  2-22  requires  the  sideband  bandwidth.  B  |.  The  derivation 

of  Equation  3-3  made  use  of  the  jammers’  fiat  spectrum  and  that  the  mainband  and  sideband  channel  spec- 
trums  are  assumed  to  have  the  same  shape  but  different  bandwidths.  As  expected,  the  jarruner’s  JNR  and 
AOA  information  is  left  unchanged  from  the  mainband,  although  the  covariance  matrix  is  modified  due  to 
the  change  in  carrier  frequency  and,  to  a  lesser  degree,  due  to  the  different  bandwidths.  If  the  mainband  and 
sideband  channel  spectrum  shapes  are  different,  then  Equation  3-3  would  remain  the  same  except  that  Ajj. 
would  require  modification  to  account  for  the  different  correlation  functions. 
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3.3  Sidcband-to-Mainband  Adaptive  Weight  Compatibility  and  Compensation 

Calculation  of  the  ASP  weights  via  the  Decoupled  Canceler  algorithm  of  Equation  2-57  using  the 
sideband  data  covariance  matrix  of  Equation  3-3  can  result  in  performance  degradation,  because  DOA 
information  is  encoded  based  on  different  carrier  frequencies  resulting  in  different  normalized  spatial  fre¬ 
quencies  for  the  jammers.  To  illustrate  this  degradation,  the  SINR  Loss  for  a  system  employing  the  side¬ 
band  data  directly  with  the  Decoupled  Canceler  algorithm  is  shown  in  Figure  3-2  and  Figure  3-3.  For  this 


example,  steering  and  beamforming  is  accomplished  using  componenu  which  have  the  same  response  in 
the  sideband  as  the  mainband.  This  is  not  an  unreasonable  assumption  since  the  beamformer  being  used  is 
the  Butler  matrix,  which  is  commonly  implemented  using  wideband  components,  and  the  steering  is 
accomplished  using  RF  phase-shifters  which  can  also  be  relatively  wideband.  A  narrow  sideband  band¬ 
width  was  chosen  (i.e.  0.1  MHz)  and  a  sideband  offset  of  4  MHz  was  used.  It  is  evident  that  the  ASP 
weights  generated  using  the  sideband  data  directly  are  incompatible  for  jammer  suppression  in  the  main- 
band.  The  problem  can  be  visualized  by  plotting  the  mainband  ASR  pattern  nulls.  The  ASR  nulls  are 
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Normalized  Target  Doppler 

Figure  3-3:  SINK  Loss  with  SIdelobe  Region  Jammers  using  Uncompensated  Decoupled  Canceler:  4 

Aux/Scenario  0 


located  at  spatial  frequency  locations  corresponding  to  the  jammer  locations  for  the  sideband  earner.  How¬ 


ever.  for  mainband  data,  these  spatial  frequencies  do  not  correspond  to  the  jammer  locations.  Therefore,  the 


resulting  nulls  in  the  mainband  are  slightly  shifted  away  from  the  jammer  locations,  and  because  the  nulls 


are  relatively  narrow,  a  significant  loss  in  jammer  suppression  can  result.  Note  that  the  loss  is  small  in  the 


case  of  the  CFB  system  which  uses  sidelobe  responses  to  cancel  the  sidelobe  jammers.  This  can  be  attrib¬ 


uted  to  the  well  matched  sidelobe  structure  of  the  Butler  matrix  resulting  in  wide-angle  januner  suppres¬ 
sion.  This  feature  can  be  shown  to  disappear  for  even  small  beamformer  channel  errors,  and  therefore 
cannot  easily  be  maintained  in  practice.  The  results  for  the  CSB  and  CSA  systems  are  more  indicative  of 
practical  performance.  In  Section  3.3.1  through  Section  3.3.4  various  compensation  techniques  are 
described,  and  performance  comparisons  for  the  case  of  sidelobe  jammers  are  provided. 


3J.1  Beamsteering  Matrix  Design 


As  a  first  step  in  the  compensation  of  the  sideband  data,  an  advantageous  beamsteering  transforma 
tion  is  selected.  Note  that  when  the  array  is  steered  off  broadside,  the  mainband  and  sideband  beams  will 
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shift  different  amounts  if  the  phase  shifters  have  the  same  response  at  both  carriers.  It  is  desirable  to  have 
the  mainbeam  peaks  at  the  different  carriers  line  up  as,  will  be  demonstrated.  Thus,  the  diagonal  beamsteer¬ 
ing  transformation  is  chosen  to  be  different  for  the  two  carriers  (i.e.  L|  ^  L^),  and  instead  the  phase  shifts 
at  the  sideband  carrier  are  designed  to  provide  the  same  steering  angle  as  the  mainband.  Since  the  sideband 
bandwidth  for  practical  systems  will  typically  be  very  narrow,  it  is  not  expected  that  this  will  cause  a  prob¬ 
lem  in  the  phase-shifter  implementation.  The  performance  improvement  for  a  system  employing  this  steer¬ 
ing  method  is  demonstrated  in  Figure  3-4.  This  initial  performance  improvement  is  compounded  in  the 


Canceler:  4-Aux/Scenario  D 

compensation  techniques  yet  to  be  discussed  since  they  can  benefit  from  the  use  of  this  beamsteering 
method. 

3  Data  Domain  Compensation  via  Focussing  Matrices 

Compensation  in  the  data  domain  refers  to  applying  a  transformation  to  the  digitized  sideband  data 
in  order  that  the  covariance  matrix  of  the  transformed  jammer  data  approximates  the  mainband  jammer 
covariance  matrix.  This  investigation  will  concentrate  on  using  a  linear  transformation,  F,  to  accomplish 
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the  desired  compensation.  Referring  to  Equation  3*3,  Equation  2-58,  and  Equation  2-24,  an  F  is  sought 
which  satisfies 


f"(L,B,)V,  (f,))  =  for  k  =  I . Nj 

or  f"(L,B,)"s,  =  (LqBo)  Sq 

where  the  columns  of  the  matrices.  Sq  and  S|,  are  the  spatial  steering  vectors  for  each  jammer  at  the  appro¬ 
priate  carrier.  For  simplicity,  the  bandwidth  decorrelation  effects,  represented  by  Aj)^  in  the  jammer  covari¬ 
ance  matrices,  have  been  ignored  in  this  equation.  It  is  easy  to  show  that  if  Equation  3-4  is  satisfied,  then 
the  jammer  covariance  matrix  for  the  compensated  sideband  data  will  be  the  same  as  for  the  mainband  jam¬ 
mer  covariance  (to  within  a  scale  factor)  if  the  bandwidth  decorrelation  effects  are  neglected.  This  proce¬ 
dure  of  using  a  linear  transformation  to  modify  data  from  one  band  in  order  to  be  used  for  processing  data 
from  a  different  band  is  closely  related  to  coherent  signal-subspace  processing  which  has  been  applied  to 
the  AOA  estimation  problem  for  w'ideband  sources.  The  resulting  linear  transformations  have  been  called 
‘focussing  matrices’  and  this  terminology  is  used  here.  Detailed  information  on  this  topic  can  be  found  in 
(231  and  (24]. 

In  this  section,  data  compensation  will  be  performed  with  the  jammer  locations,  power,  and  num¬ 
ber  being  unknown.  Since  a  jammer  can  presumably  be  located  at  any  angle,  the  focussing  matrix  is 
designed  for  a  set  of  prespecified  angles  evenly  distributed  in  spatial  frequency  (of  the  mainband)  over  the 
entire  field  of  view.  Thus  the  Sq  matrix  will  not  correspond  to  the  unknown  jammer  locations,  but  will 
instead  reflect  the  preselected  locations  of  calibration  sources.  Accordingly,  the  S|  matrix  corresponds  to 
the  same  source  angles,  but  will  have  a  structure  reflecting  the  sideband  carrier.  In  this  way,  the  focussing 
matrices  are  precalculated  using  receiver  calibration  data.  The  measurements  would  allow  the  matrix  prod¬ 
uct  (LB)^D  to  be  determined  for  the  mainband  and  sideband  receivers.  The  approach  used  here  is  to  use  as 
many  focussing  angles.  N,.  as  possible  while  allowing  Equation  3-4  to  have  an  exact  solution.  Note  that  the 
matrices  Si  (N^xN,).  L,  (N^-xN^).  and  B,  (N^xN^)  are  all  full  tank  for  any  practical  system  design.  If  we 
require  N,  ^f,  then  it  is  straight  forward  to  show  that  an  exact  solution  to  Equation  3-4  is 
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f”  =  (LoBo)”SoI  (L,B,)”s,l  .  (-V5) 

where  the  superscript  *+’  denotes  Pseudo- In  verse  (i.e.  Moore-Penrose  inverse).  The  solution  of  Equation  3- 
5  is  exact  whenever  the  term  in  *[]’  brackets  is  full  column  rank.  For  the  results  of  this  section,  the  number 
of  focussing  angles,  N,,  and  focussing  channels,  Nf.  are  chosen  to  be  equal  to  the  number  of  sensor  chan¬ 
nels  (i.e.  N4=N(=Nc).  If  it  is  desired  to  focus  at  more  angles  than  Nf  (i.e.  N,  >Nf),  then  it  can  be  shown  that 
the  solution  given  in  Equation  3-5  is  also  the  solution  to  the  more  general  problem  of  finding  the  matnx.  F. 
which  minimizes  the  Frobenius  norm  of  the  error  matrix  (i.e.  the  matrix  difference  of  the  two  sides  of 

Equation  3-4). 

Although  the  matrices  Bq,  Lq-  L,  have  already  been  specified,  there  is  some  room  for  selecting 
the  sideband  beamformers.  B,.  for  each  system.  In  order  to  keep  the  effects  of  sideband-to-mainband  chan¬ 
nel  mismatch  to  a  minimum,  it  is  desirable  to  have  the  beamformers  matched.  Also  to  keep  system  co.st 
down,  the  sideband  beamformer  should  use  as  much  of  the  mainband  hardware  as  possible.  Hence,  for  the 
CFB-JDL  and  the  CSB-JDL  systems  the  sideband  beamformer  is  chosen  to  be  a  full  Butler  matrix,  whereas 
for  the  eSA-JDL  the  same  auxiliary  channels  are  used  as  well  as  the  JDL  beams,  plus  enough  additional 
Butler  beams  (adjacent  to  the  JDL  beams)  in  order  to  attain  N^  sideband  channels.  Performance  results  for 
this  technique  using  a  fixed  set  of  calibration  source  locations  is  given  in  Figure  3-5.  Note  that  although  the 
CSB-JDL  and  CS  A-JDL  systems  show  improvement,  the  Cre-JDL  suffers  some  degradation  in  the  central 
doppler  region.  The  degradation  can  be  attributed  to  the  calibration  source  locations  not  being  aligned  with 
the  Butler  beam  peaks  when  the  array  is  steered  off  broadside  and  not  aligned  with  a  calibration  angle. 
Another  troublesome  feature  of  this  technique  is  that  the  focussing  matrices  are  functions  of  the  array  steer¬ 
ing  angle,  and  therefore  a  set  of  focussing  matrices  are  needed  to  cover  a  range  of  steering  angles.  Both  of 
these  undesirable  attributes  can  be  abrogated  by  modifying  the  calibration  method  as  explained  below. 

As  an  alternative  to  the  Focussing  calibration  discussed  above,  consider  the  following.  Instead  of 
requiring  a  fixed  set  of  calibration  angles,  the  angle  locations  are  shifted  in  unison  as  the  array  steenng 
angle  changes  such  that  the  central  source  is  kept  aligned  with  the  central  target  beam  (i.e.  mainbeam).  For 
broadside  steering  the  sources  are  chosen  as  before,  which  are  aligned  with  the  peaks  of  the  Butler  beams 
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Figure  3-5:  SINR  Loss  with  Sidelobe  Region  Jammers  using  Steering-Dependent  /  Focussing- 
Compensated  Decoupled  Canceler:  4-Aux/Scenario  D 

with  N4=N<.=Nf.  Now  making  use  of  the  beamsteering  definition  of  Section  3.3.1,  the  fixed  steering  matri¬ 
ces  of  Equation  3-4,  Sq  and  Sj,  are  replaced  by  steering  matrices  which  are  a  function  of  the  steering  angle. 


I>0  =  LqSq 


D,  =  L,S, 


and  which  keep  the  mainband  and  sideband  mainbeams  aligned.  Substituting  these  steering  matrices,  Dj 
and  D|,  for  Sq  and  Sjin  Equation  3-4,  and  using  the  unitary  feature  of  the  beamsteering  matrix,  the  focus¬ 
sing  matrix  equation  to  be  solved  becomes 


=  B„"s„. 


v.1iich  is  now  independent  of  the  array  steering  angle  and  reflects  a  broadside  steering  condition.  Clearly, 
the  solution  to  Equation  3-7  has  the  same  form  as  Equation  3-5,  but  with  the  beamsteering  matrices  absent. 
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The  implication  is  that  as  long  as  the  beamsteering  definition  of  Section  3.3. 1  can  be  implemented  accu¬ 
rately.  only  a  single  focussing  matrix  needs  to  be  calculated  using  only  the  fixed  broadside  set  of  calibration 
sources.  Performance  results  given  in  Figure  3-6  using  this  method  for  calculation  of  the  focussing  matrix 


Normalized  Target  Doppler 

Figure  3-6:  SINR  Loss  with  Sidelobe  Region  Jammers  using  Steering-Independent  /  Focussing- 
Compensated  Decoupled  Canceler:  4-Aux/Scenario  D 

shows  an  improved  performance  for  the  CFB-JDL  alluded  to  above.  The  CSB-JDL  system  also  shows 
slight  improvement,  whereas  the  CSA-JDL  system  shows  minor  degradation  at  some  dopplers.  Focussing 
at  the  Butler  peaks  evidently  is  beneficial  to  systems  employing  butler  beams  for  jammer  suppression, 
while  a  system  using  omni-directional  auxiliaries  does  not  appear  to  benefit. 

3  Spectral  Analysis  Aided  Focussing  Compensation 

It  is  of  interest  to  determine  how  much  the  data  domain  focussing  method  can  be  improved  if  the 
focussing  matrix  is  calculated  adaptively,  after  spatial  spectral  analysis  has  located  the  jammers.  There  are 
many  good  techniques  available  in  the  literature  for  AOA  processing  of  narrowband  data  (e.g.  see  (25)  for 
element  space  and  [26]  for  beamspace  techniques),  and  it  is  assumed  here  that  perfect  estimates  are  avail¬ 
able.  The  focussing  matrix  is  designed  using  the  steering-independent  method  and  focusing  at  the  peaks  of 


the  cluiter/target  beams  and  at  the  jammer  locations.  Thus,  the  Sq  used  in  Equation  3-7  is  generated  using  a 
steering  matrix.  Dq.  which  has  Nj  columns  corresponding  to  the  clutter  beams,  and  additional  columns 
corresponding  to  the  jammer  locations  transforming  according  to  Equation  3-6.  As  expected,  the  perfor¬ 
mance  results  are  quite  good  as  shown  in  Figure  3-7  and  Figure  3-8.  It  should  be  noted  that  in  the  literature. 


Normalized  Target  Doppler 

Figure  3-7:  SINR  Loss  with  Sidelobe  Region  Jammers  using  Spectrally-Assisted  Focussing- 
Compensated  Decoupled  Canceler:  4-Aux/Scenario  D 

when  a  transformation  precedes  the  AOA  algorithm  it  often  is  required  to  be  unitary  in  order  to  leave  the 
noise  subspace  unperturbed.  This  is  not  a  problem  for  the  CFB-JDL  and  CSB-JDL,  but  an  alternative  side¬ 
band  beamformer  (e.g.  the  Identity  implying  no  beamforming)  may  be  needed  for  the  CS A-JDL  system. 

Some  apparent  drawbacks  of  AOA  assisted  focussing  are  ( I)  significant  amount  of  additional  com¬ 
putation  for  the  AOA  prbeessing,  and  (2)  many  more  calibration  points  are  required  to  achieve  full  poten¬ 
tial.  However,  once  the  AOA  information  is  obtained,  new  opportunities  for  performance  improvement  arc 
available.  For  instance,  AOA  information  can  identify  clutter  and  signal  region  jammers,  thereby  allowing 
dynamic  allocation  of  spatial  DOF  for  those  interferers  (e.g.  force  the  ASP  to  ignore  them  and  increase  the 
spatial  DOF  out  of  the  ASP  in  order  to  process  these  jammers  in  the  STAP). 
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Figure  3-8:  CSB-JDL  ASR  Pattern  y>ith  Sidelobe  Region  Jammers  using  Spectrally-Assisted 
Focussing-Compensated  Decoupled  Canceler:  4-Aux/Scenario  D 


33.4  Data  Domain  Compensation  via  Dual-Sideband  Combining 

The  next  form  of  compensation  to  be  discussed  attempts  to  shift  and/or  widen  the  ASP  null  by 
combining  two  symmetrically  located  sidebands  (with  respect  to  the  mainband)  prior  to  weight  calculation. 
The  approach  used  here  is  to  have  both  sideband  channels  designed  the  same  way  with  regard  to  their 
respective  carrier  frequencies  and  to  minimize  the  duplication  of  hardware.  The  beamsteering  unit  is 
designed  as  in  Section  33.1,  but  now  having  an  additional  phase  shift  specification  for  the  added  sideband. 
The  beamformers  are  assumed  wideband  enough  to  have  the  same  response  for  both  sidebands.  A  simple 
modification  to  the  SBC  block  diagram  of  Figure  3-1  can  be  used  to  obtain  the  combined  sideband  WNJ 
signals  using  a  single  sideband  filler /down-conversion /A/D  chain.  Prior  to  the  sideband  filler,  the  signals 
are  down-converted  so  that  the  mainband  is  at  baseband.  Note,  this  is  tantamount  to  using  the  mainband 
data  just  prior  to  the  Mainband  Filter  block,  because  the  beamformers  arc  designed  to  be  the  same.  Then  a 
passband  sideband  filter  is  applied,  followed  by  down-conversion  using  the  sideband  carrier  offset  fre¬ 
quency,  f^ff,  and  finally  A/D  conversion.  The  beamformers  do  not  need  to  be  the  same  to  use  this  technique. 
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however,  more  hardware  will  be  required,  and  the  system  w'ill  be  more  susceptible  to  mismatch  between 
the  sideband  and  mainband  jammer  data. 

The  dual-sideband  data  covariance  matrix  that  is  being  generated  is  seen  to  be  the  summation  of 
two  covariance  matrices  which  have  the  same  form  given  in  Equation  3-3,  but  corresponding  to  different 
carrier  frequencies.  In  general  then,  it  is  assumed  the  data  from  each  sideband  is  uncorrelated,  and  the 
resulting  covariance  matrix  is 


=  (L,B,) 


(LjBj) 


(L,B,)  -H 

(LjBj) 


(3-8) 


However  for  the  specific  systems  considered  here,  by  design  the  sideband  bandwidths,  beamformers  and 
carrier  frequencies  are  given  by 


B,  =  Bj. 


B,=B2. 


^  ~  ^off 


and 


(3-9) 


Performance  results  after  applying  this  dual  sideband  model  to  each  of  the  baseline  systems  is  given  in  Fig¬ 
ure  3-9.  Comparing  these  results  with  those  of  Figure  3-7  show  that  the  Spectrally-Assisted  Focussing 
Matrix  produces  slightly  better  SINR  performance  for  the  CSB-JDL  and  CSA-JDL  systems  even  though 
the  Dual-Sideband  approach  provides  better  jammer  suppression  as  indicated  by  the  JNR  levels  given  in 
the  plot  legends.  This  difference  is  pointed  out  to  illustrate  that  in  addition  to  the  jammer  suppression. 
STAP  performance  is  also  dictated  by  the  clutter  subspace  at  the  output  of  the  ASP.  The  Dual-Sideband 
approach  causes  more  perturbation  to  the  clutter  subspace  in  an  effort  to  generate  a  wide  null  as  can  be  seen 
by  comparing  the  ASR  patterns  in  Figure  3-10  and  Figure  3-8  for  example.  Note  that  the  clutter  subspace 
width  has  been  decreased.  For  systems  where  the  sidebands  are  close  together  this  may  not  be  a  problem. 
However,  as  the  sidebands  move  apart  the  clutter  subspace  will  continue  to  degrade.  In  this  case  a  tech- 
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Decoupled  Canceler:  4-Aux/Scenario  D 


nique  such  as  Space-Range  processing  (comparable  to  the  so-called  Frost  array  which  uses  tapped-dclay 
lines,  see  e.g.  [27])  may  be  used  in  order  to  generate  a  wide  null  without  degrading  the  clutter  subspace. 
This  topic  is  left  for  future  investigation. 

3.3.5  Weight  Domain  Compensation  via  Double-Sideband  Weight  Averaging 

In  this  section  a  different  compensation  technique  utilizing  two  sidebands,  called  Double-Sideband 
Weight  Average  Compensation  (DWAC),  is  discussed.  The  fundamental  concept  is  to  generate  ASP  weight 
vectors  separately  for  each  sideband  and  then  combine  them  to  produce  a  weight  vector  appropriate  for 
jammer  cancellation  in  the  mainband.  Digitized  data  sets  for  two  sidebands  located  symmetrically  about 
the  radar  mainband  are  generated.  ASP  weights  for  each  sideband  produce  nulls  on  opposing  sides  of  the 
desired  mainband  location.  The  intent  is  to  combine  them  in  a  way  such  that  the  spatial  responses  combine 
to  place  the  null  at  the  desired  location  as  in  the  Dual-Sideband  compensation  method  of  Section  3.3.4,  but 
without  the  additional  bandwidth  effects  inherent  in  that  approach.  Of  course,  the  data  sets  could  be  gener¬ 
ated  by  using  two  separate  sideband  receivers  constructed  as  in  Figure  3-1,  but  more  efficient  ways  which 
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Figure  3-10:  CSB-JDLASR  Pattern  with  Sidelobe  Region  Jammers  using  Dual-Sideband 
Compensated  Decoupled  Canceler:  d-Aux/Scenario  D 

use  the  same  hardware  for  both  channels  can  easily  be  developed.  For  instance,  the  down-conversion  stage 
before  the  Sideband  Filter  can  alternately  use  two  opposing  sideband  offset  frequencies,  If^ff,  thereby  gen¬ 
erating  digitized  data  from  alternating  sidebands.  Thus  the  DWAC  method  averages  the  weight  vectors 
generated  from  the  two  sideband  data  sets  using  the  same  cancelation  algorithms.  Thus  in  the  case  of  the 
Decoupled  Canceler  of  Equation  2-57,  sample  covariance  matrix  estimates  (corresponding  to  Equation  3- 
3)  for  each  data  set  are  calculated  and  the  DWAC  compensated  weight  matrix  is  given  by 


W  =  i-  (Wj+Wj)  . 


(3-10) 


Note  that  the  gain  constraint  of  the  DC  is  preserved  since  both  W|  and  W2  satisfy  the  same  constraint.  In 
general,  one  should  verify  that  the  weight  vector  averaging  makes  sense  for  the  cancelation  algorithm  being 
used  in  terms  of  satisfying  any  constraints  and  interference  suppression.  For  the  cancelation  algorithms 
considered  in  this  investigation,  it  is  found  that  the  desired  null-shift  is  indeed  obtained  by  this  averaging 
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technique.  Cleaxly.  the  scale  factor  of  1/2  is  superfluous,  but  left  here  to  illustrate  the  weight  balancing 
being  prop>osed.  Very  acceptable  performajice  is  obtained  from  this  simple  technique  as  evidenced  by  the 
JNR  and  SIKR  Loss  results  given  in  and  Figure  3- 1 1 .  It  is  also  instructive  to  look  at  the  JNR  residual  as  .n 


Figure  3-11:  SINR  Loss  for  Sidelobe  Region  Jammers  using  Double-Sideband  Average  Compensated 

Decoupled  Canceler:  4-Aux/Scenario  D 

function  of  sideband  offset  frequency  given  in  Figure  3-12.  The  Double-Sideband  Averaging  provides 
fairly  uniform  performance  over  a  ±5  MHz  offset  range.  Although  not  shown  here,  the  Focussing  mairi.v 
approaches  (with  the  exception  of  the  spectrally-assisted  version)  do  not  exhibit  this  uniform  JNR  perfor¬ 
mance.  but  rather  have  a  minimum  at  zero  offset  frequency,  and  can  degrade  rather  quickly  away  from  zero 
offset.  In  a  sense  then,  the  double-sideband  method  can  be  regarded  as  having  larger  “compensation  band¬ 
width”. 

3  J.6  Summary  of  Sideband  Compensation  Techniques 

In  this  section,  a  sideband  method  was  described  for  acquisition  of  jammer-only  secondary  data, 
along  with  several  methods  for  using  the  sideband  data  to  adaptively  generate  spatial  weights  for  the  sup¬ 
pression  of  mainband  jamming.  The  compensation  techniques  considered  belong  to  two  general  categories: 
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Decoupled  Canceler:  4-Aux/Scenario  D 

(I)  data  domain  compensation  whereby  the  sideband  data  is  transformed  to  take  on  the  characteristics 
indicative  of  the  mainband  interference  prior  to  weight  calculation,  and  (2)  weight  domain  compensation 
where  spatial  weight  vectors  generated  from  uncompensated  sideband  data  are  modified  to  provide  the 
desired  interference  cancelation  in  the  mainband.  Among  the  data  domain  techniques,  the  spectrally- 
assisted  focussing  matrix  approach  shows  the  best  performance  potential.  However,  because  of  it’s  addi¬ 
tional  requirement  for  pre-location  of  the  jammers,  the  steering  independent  focussing  matrix  approach 
with  minor  additional  SINR  Loss  may  be  preferred.  The  Dual-Sideband  approach  also  shows  acceptable 
performance,  but  in  systems  employing  two  sidebands  the  Double-Sideband  Averaging  method  may  be 
preferred  since  it  exhibits  slightly  better  performance  and  does  not  consume  spatial  DOF  to  generate  wide 
spatial  nulls.  However,  it  should  be  noted  that  wide  spatial  nulls  can  provide  protection  against  mainband- 
to-sideband  channel  mismatch,  but  this  presupposes  the  use  of  space-range  processing,  which  is  beyond  the 
scope  of  this  report. 
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4  Summary 


4.1  Conclusions 

Several  cascaded  archilectures  have  been  proposed  as  an  implemenlation  of  (he  SAPS  concept. 
The  performance  results  for  these  cascaded  systems  compare  favorably  with  the  simultaneous  STAP 
approaches,  and  verify  their  viability  in  future  radar  designs.  A  key  component  of  the  WNJ  pre-suppression 
concept  is  the  acquisition  of  jammer-only  secondary  data.  The  sideband  technique  described  herein  has  the 
unique  advantage  over  other  techniques  that  it  allows  collection  of  secondary  data  during  any  portion  of  an 
RPI.  Sideband  compensation  methods  were  proposed  and  analyzed  showing  that  ample  Jammer  suppres¬ 
sion  can  be  achieved  in  the  mainband  while  using  the  sideband  secondary  data  for  weight  calculation. 

4.2  Further  Work 

The  work  done  so  far  in  the  area  of  jammer  pre-suppression  in  conjunction  with  STAP  has  only 
scratched  the  surface.  More  investigation  into  dealing  with  jammers  in  the  clutter  and  target  regions  a.s  well 
as  other  sideband  compensation  methods  is  prudent.  Also,  studying  the  effects  of  imperfect  beamformers 
(i.e.  channel  errors)  and  mainband-lo-sideband  channel  mismatch  is  of  great  practical  interest.  A  well 
defined  wholistic  approach  to  the  design  of  “optimar  beamfoimer  transformations  is  another  possible  area 
to  improve  the  proposed  systems.  Perhaps  more  importantly,  is  the  expansion  of  (he  divide-and-conquer 
concept  for  handling  interference  generated  by  sources  with  dissimilar  characteristics.  When  faced  with 
real-world  problems  which  can  often  present  multiple  forms  of  interference,  being  able  to  focus  ones 
resources  on  each  separately  can  be  very  advantageous.  One  example,  for  which  the  sideband  concept 
described  here  may  provide  a  skillful  solution,  is  the  problem  of  jammer  multipath  (also  known  as  **hor 
clutter).  Sideband  data  could  be  used  to  supply  a  Space-Range  Processor  (SRP)  with  a  sufficient  amount  of 
secondary  data  to  allow  the  application  of  a  sizeable  adaptive  space-range  algorithm. 
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Appendices 


Appendix  A:  Optimum  Receiver  for  Random  Target  Amplitude  Model 

In  this  appendix,  a  sufficient  statistic  for  the  binary  hypothesis  test  of  Section  2.2.1  is  generated  by 
algebraic  manipulation  of  the  likelihood  ratio.  The  resultant  statistic  is  the  optimal  Neyman-Pearson 
receiver  under  the  assumed  Gaussian  data  distribution  models.  The  data  model  for  the  pxl  vector,  x.  is 
given  by 

X  =  as  +  n  where  n  ~  CN(0,  R)  and  under  each  hypothesis 

2  ’  ^  ^ 
=  0  and  H|:a~CN(0,a) 

where  is  read  “is  distributed  as”  and  “CN”  represents  the  complex  normal  distribution  with  probability 
density  function  (see  e.g.  [191) 

pfxjm,  R)  =  Jt'^’lRl  'exp[-(x-in)^R  '  (x-m)  1  .  (A2) 

It  is  comp)oscd  of  a  target  vector  with  random  complex  amplitude,  a,  and  specified  structure,  s,  embedded 
in  additive  zero-mean  (circularly)  complex  Gaussian  interference.  It  is  assumed  that  the  interference,  n,  and 
the  target  amplitude  are  independently  distributed.  The  data  has  zero-mean  under  both  hypotheses,  but  the 
covariance  matrix  is  different  with 

Hq:E{xx^}  =R  and  Hj:E{xx  }  =R  +  o  ss  =Rj.  (A3) 

Under  hypothesis  Hq,  the  data  obviously  has  the  same  distribution  as  the  interference.  Under  Hj  it  is  easy 
to  show  that  the  data  is  also  CN  distributed  as  follows.  Since  the  interference  is  CN.  any  linear  combination 
is  also  CN  (see  [19]  p.  375).  Therefore,  by  the  independence  of  a  and  n,  for  any  pxl  vector,  a,  by  Proposi¬ 
tion  9.15  of  [19] 
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a”x  -CN(0.a“R,a). 


Thus  by  Definition  9.2  of  [19]  the  data  is  distributed  as 


H,:x-CN(0,R,) 


Now  the  likelihood  ratio  statistic  is  given  by 


PxlH,  iRlexp  l-x^R, 'xj 
PxIHq  |R||exp[-x^R  'xj 


By  taking  the  logarithm  of  /  and  ignoring  terms  which  do  not  depend  on  the  data,  an  equivalent  statistic  is 


Ti'  =  x”(R“'-R;')x. 


Using  the  matrix  inversion  !emma  (a  general  form  can  be  found  in  [20])  and  Equation  A3  it  is  found  that 


p-l  _  „-i  R~*ss”r  * 

2  H„-l 
o  +  s  R  s 


from  which  the  statistic  can  be  written  as 


.  x”r"'ss”r"'x 

^  =  2  h^-T- 

o  +s  R  s 


Now  the  likelihood  ratio  decision  rule  can  be  written  as 


I  p  ^1' 

Ti  =  Is  R  xl  ^  TIq 

HO 


(AID) 
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which  is  seen  to  involve  colored-noise  matched  filtering  of  the  input  data. 

Appendix  B:  Simplification  of  the  GLR  Conditional  Detection  Probability 

In  this  appendix,  the  conditional  detection  probability  formula  for  the  GLR  detector  given  by 
Equation  2-37  is  simplified.  The  formula  for  P^ip  is  rev/ritten  here  as 


j  =  1  m  =  C 


m 


=  0  Iq  ^  I 


m  +  I 


where  the  following  substitutions  have  been  used; 


(Bl) 


R  -  PPjC*  ~^oT^ 


.and  M  =  K-N,-t-  I  . 

•  -Hot 


(B2) 


The  second  summation  in  Equation  B I  can  be  put  in  the  form  of  a  geometric  progression.  Then  using  the 


closed  form  expression  given  in  [21]  p.  1  yields 


j  -  ^  m 

X— ^ 

m  =  0  Iq  +  '  1 


m  +  I 


(B3) 


Substituting  the  result  of  Equation  B3  into  Equation  B 1  it  is  seen  that  Pjip  can  be  written  as  the  sum  of  two 
finite  binomial  power  series.  Using  the  closed  form  expression  for  a  binomial  power  series  given  by  [2 1  ]  p. 
21  produces  the  formula 


j  =  1  i  =  o 

Substituting  Equation  B3  and  Equation  B4  into  Equation  Bl  results  in  the  following  simplification  for  P^ip 
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^p-'-Pr 

j  =  I 

=  l-P,  (B5) 

-  p .  r__ppTi!_r 
^  LpPt(‘-^ot)  +  U 


Now,  using  the  relationship  between  the  threshold  and  Pf  given  by  Equation  2-35  produces  the  desired  final 
form  for  the  conditional  detection  probability 


pPt+  I 


pPx  +  Pf 


■I/M 


M 


(B6) 


Appendix  C:  Closed  Form  Expression  for  the  Integrated  System  Response 

A  closed  form  expression  for  the  integrated  System  Response,  Aj- defined  by  Equation  2-54  is 
derived  for  the  case  of  data  uniformly  sampled  in  space  and  time.  Under  these  conditions,  the  expressions 
for  the  spatial  and  temporal  steering  vectors,  s,  and  Sj,  of  Equation  3- 1 1  are  appropriate.  Using  Equation  2- 
49.  Equation  2-54  can  be  written  as 


fts  I  ss  s  \ 


6^M,  f„)  = 


f  -Af/2f  -Af/2 
1$  r  st  t 


H 

^of^T 


(Cl) 


where  the  terms  which  are  not  a  function  of  the  spatial  or  Doppler  frequencies  can  be  taken  outside  the  inte 
gral.  Furthermore,  using  the  properties  of  the  Kronecker  product  and  separating  sole  ftinctions  of  spatial 
and  Doppler  frequencies  yields  the  following  formula 
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(C2) 


Notice  that  each  integral  is  of  the  same  form  with  differences  only  in  the  parameters  of  the  integrand  and 
the  limits  of  integration.  Also  note  that  the  integrands  are  Hermiiian-Toeplitz  matrices  so  that  only  the  ele¬ 
ments  of  the  first  column  need  to  be  integrated.  For  convenience,  rewrite  Equation  C2  as 

AT(Af,.  -  (C3) 

where  the  corresponding  definitions  for  the  matrix  integrals,  Aj  and  Aj,  are  obvious.  The  elements  of  A,  and 
Aj  can  be  evaluated  using  elementary  integration  formulas.  The  resulting  expressions  are 

Aj  =  Toeplitz(Sjj  •  d^)  2md  A^  =  Toeplitz(Sj|  •  d^) 
where  the  elements  of  d^  are  d^fn)  =  (Af^)  sine  (nAf^)  ,  n  =  0,  ....  -  1  .  (C4) 

and  similarly  for  d^  ,  d^fn)  =  (Af^)  sine  (nAf^)  ,  n  =  0, ....  -  1 

The  function  ToeplitzO  represents  a  Hermitian-Toeplitz  matrix  with  the  first  column  given  by  it's  argument, 
and  the  sinc()  function  is  defined  in  the  usual  manner  as 

sine  (x)  s  sin  (tix) /itx .  (C5) 

The  heavy  dot,  represents  a  matrix  Hadamard  product  (i.e.  element-by-element  product). 
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